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Abstract 

The  mechanical  response  of  PMR-15  neat  resin  was  investigated  at  288  °C.  The 
effect  of  loading  rate  on  monotonic  stress-strain  behavior  was  explored  in  monotonic 
tests  at  several  constant  stress  rates.  Considerable  rate  dependence  was  observed, 
especially  on  the  unloading  path.  Effect  of  prior  stress  rate  on  creep  behavior  was 
evaluated  in  creep  tests  preceded  by  uninterrupted  loading  to  a  target  stress  where  loading 
rate  was  changed  from  test  to  test.  Creep  response  was  dependent  on  the  prior  stress  rate. 
Effect  of  loading  history  was  studied  in  stepwise  creep  tests,  where  specimens  were 
subjected  to  a  constant  stress  rate  loading  followed  by  unloading  to  zero  stress  with 
intermittent  creep  periods  during  both  loading  and  unloading.  Comparison  of  creep 
strains  accumulated  during  a  stepwise  creep  test  to  those  accumulated  during  creep 
preceded  by  uninterrupted  loading  indicate  that  the  prior  stress  history  affects  the  creep 
behavior.  Negative  creep  as  well  as  creep  rate  reversal  was  observed  on  the  unloading 
path. 

A  nonlinear  viscoelastic  model  (Schapery’s  fonnulation)  was  characterized  using 
creep  and  recovery  tests  at  288  °C.  The  model  was  verified  by  comparing  the  predictions 
with  experimental  results  obtained  in  monotonic  loading/unloading  tests  and  single-step 
as  well  as  multi-step  creep  tests.  The  model  qualitatively  predicted  creep  response  to 
single-  and  multi-step  creep  tests,  including  negative  creep  and  creep  rate  reversal  during 
unloading.  However,  predictions  were  not  quantitatively  accurate.  The  model  was 
unable  to  accurately  predict  the  recovery  behavior  and  could  not  account  for  rate  effects. 
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SOME  ASPECTS  OF  THE  MECHANICAL  RESPONSE  OF  PMR-15  NEAT 


RESIN  AT  288°C:  EXPERIMENT  AND  MODELING 


I.  Introduction 

The  aerospace  industry’s  unrelenting  drive  for  the  advancement  of  aircraft 
performance  simultaneously  fuels  the  need  for  high  performance  structural  materials. 
Composite  materials  stand  at  the  forefront  of  this  search  as  their  high  strength  and 
stiffness  to  weight  ratios  offer  numerous  advantages  over  other  materials.  Composites 
materials  are  not  a  novel  concept.  Modern  composite  materials  first  played  a  substantial 
role  in  military  applications  with  the  inclusion  of  boron-reinforced  epoxy  composites  in 
the  skins  of  the  U.S.  Navy’s  F-14.  While  initially  composite  materials  were  only  used  in 
secondary  structures,  the  evolution  of  aircraft  design  has  already  witnessed  substantial 
advancements  with  the  F-22  Raptor  carrying  24%  of  its  weight  in  thermoset  composites. 

Polymer  matrix  composites  (PMCs)  have  been  incorporated  into  aircraft  structural 
components  for  many  decades.  In  the  1980’s,  many  low  speed  aircraft  such  as  the 
Boeing  757  and  McDonnell  Douglas  MD-80  series  contained  components  made  of  epoxy 
resins  reinforced  with  aramid,  carbon,  and  E-glass  fibers  (9:36).  While  these  composites 
met  the  requirements  of  the  aforementioned  aircraft,  the  epoxy  resins  were  limited  to 
operating  temperatures  below  130°C.  Therefore,  researchers  sought  out  more  aggressive 
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polymer  systems  that  could  withstand  temperatures  beyond  that  of  the  epoxy  resins.  The 
result  was  intense  investigation  of  a  group  of  polymers  called  polyimides  which  are 
recognized  today  as  extremely  stable  and  environmentally  resistant  polymer  systems 
(30:477).  Polyimides  are  appealing  for  their  use  as  matrix  resins  for  high-temperature 
composite  applications  due  to  their  high  toughness,  high  tensile  strength,  high  modulus, 
and  resistance  to  temperature  and  solvents  (30:477).  Current  efforts  focus  on 
incorporating  lightweight  High  Temperature  Polymer  Matrix  Composites  (HTPMCs)  into 
engine  or  body  parts  of  supersonic  aircraft. 

Prior  to  their  use  in  advanced  aerospace  applications,  however,  the  behavior  of  these 
advanced  materials  under  a  variety  of  loading  conditions  must  be  investigated  and  a 
predictive  model  must  be  developed  and  validated.  Furthermore,  the  use  of  a  composite 
material  requires  the  ability  to  predict  not  only  the  composite’s  behavior,  but  the  behavior 
of  its  constituents.  For  composites  utilized  primarily  at  high  temperatures,  it  is  therefore 
necessary  to  develop  a  predictive  model  for  the  behavior  of  the  matrix  resin  at  such 
temperatures  in  effort  to  predict  the  overall  behavior  of  the  composite  itself.  The  present 
effort  investigates  the  mechanical  response  of  PMR-15,  a  polyimide  neat  resin.  This 
research  will  focus  on  the  use  of  an  existing  model  to  determine  if  it  can  be  used  to 
predict  the  behavior  of  the  polymer  neat  resin.  This  task  is  a  critical  step  in  advancing  the 
predictive  model  framework  and  use  of  HTPMCs. 

The  following  sections  will  provide  a  detailed  explanation  of  this  research.  First,  a 
background  of  composite  materials  with  a  focus  on  polymers  will  be  presented.  This 
section  will  also  provide  a  thorough  discussion  of  the  expected  mechanical  behavior  of 
polymer  materials  in  addition  to  the  modeling  efforts  that  have  been  developed  for 
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predicting  this  sort  of  behavior.  Next,  an  explanation  of  the  material  being  tested  and  the 
experimental  setup  and  procedures  used  for  testing  will  be  described.  The  presentation  of 
experimental  results  and  model  capabilities  will  follow.  Finally,  recommendations  for 
future  efforts  and  conclusions  will  be  provided. 
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II.  Background 


This  section  provides  a  general  overview  of  polymer  matrix  composites  (PMCs) 
including  a  description  of  polymers  and  their  specific  use  as  neat  resins  in  composite 
materials.  The  use  of  PMCs  in  various  aerospace  applications  will  also  be  presented. 
Next,  the  basic  concepts  of  viscoelastic  theory  used  to  describe  the  behavior  of  polymers 
will  be  reviewed.  The  corresponding  constitutive  equations  used  to  describe  materials 
which  undergo  viscoelastic  behavior  are  presented,  including  both  the  linear  and  non¬ 
linear  fonnulations.  Finally,  a  summary  of  previous  efforts  involving  PMR-15  and  the 
modeling  of  high- temperature  polymer  deformations  will  be  given  and  the  objective  of 
this  research  will  be  stated. 

Polymer  Matrix  Composites 

A  composite,  by  definition,  is  a  “material  system  consisting  of  two  or  more  phases 
on  a  macroscopic  scale,  whose  mechanical  performance  and  properties  are  designed  to  be 
superior  to  those  of  the  constituent  materials  acting  independently”  (5:1).  In  polymer 
matrix  composites  (PMCs),  one  phase  of  the  composite  is  a  polymeric  material.  PMCs 
are  generally  used  for  structural  purposes  in  engineering  applications.  They  may  assume 
a  variety  of  forms  depending  on  the  reinforcement  phase  distributed  throughout  the 
matrix.  Continuous-fiber  reinforced,  particulate-reinforced,  and  even  nanoparticle  filled 
materials  may  be  classified  as  PMCs  (3:4).  In  high-perfonnance  structural  applications, 
PMCs  typically  consist  of  micron-diameter  fibers  bound  in  a  polymer  resin  material.  The 
fibers  are  responsible  for  the  stiffness  and  strength  properties  of  the  composites  whereas 


4 


the  matrix  provides  protection  and  support  for  the  fibers.  In  lower  perfonnance 
composites,  it  is  possible  for  the  matrix  to  be  the  main  load  bearing  constituent  for  the 
material.  It  is  necessary  to  gain  an  understanding  of  the  polymers  used  as  matrices  in 
PMCs  before  undertaking  a  review  of  their  role  in  aerospace  applications. 

Polymers 

The  word  polymer  literally  means  “many  parts”.  In  fact,  a  polymer  is  a 
“macromolecule  that  contains  many  groups  of  atoms,  called  monomeric  units”  covalently 
bonded  together  (21 : 1).  Polymers  are  usually  classified  into  three  groups  based  upon 
their  thermomechanical  properties.  These  groups  are  as  follows:  thermoplastics, 
elastomers,  and  thermosets. 

Thermoplastics,  often  simply  referred  to  as  plastics,  become  fluid  when  heated  and 
can  be  easily  processed  using  techniques  such  as  injection  molding  and  extrusion  (21:6, 
14:10).  Their  structure  can  be  either  linear  or  branched.  Thennoplastics  can  be  further 
divided  into  two  classifications:  amorphous  and  crystalline.  These  divisions  are  based 
upon  the  manner  in  which  the  plastic  is  arranged  upon  cooling.  Most  thermoplastics  do 
not  crystallize  upon  cooling  as  this  requires  reorganization  of  their  highly  coiled 
macromolecules  which  become  entangled  in  the  liquid  state  (14:10).  Those  that  are 
capable  of  some  crystallization  usually  have  both  organized  and  amorphous  regions. 
These  crystalline  polymers  are  characterized  by  their  melting  temperature,  Tm.  Most 
plastics,  however,  are  incapable  of  crystallization  upon  cooling  and  assume  a  random 
structure  even  upon  annealing.  These  amorphous  polymers  are  characterized  by  their 
glass  transition  temperature,  Tg,  at  which  they  transform  suddenly  from  a  “glassy”  hard 
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state  to  a  “rubbery”  soft  state.  Below  their  Tg,  amorphous  polymers  are  “frozen”  in 
position,  but  beyond  this  temperature  limit  they  can  move  freely  (14:10). 

Elastomers,  or  rubbery  polymers  as  they  are  commonly  referred  to,  are  cross-linked 
polymers  that  can  be  easily  defonned  many  times  beyond  their  original  length.  However, 
these  polymers  rapidly  recover  their  original  dimensions  upon  release  of  the  applied 
force.  Examples  of  common  elastomers  include  rubber  and  neoprene. 

Thermoset  polymers  are  the  most  predominant  types  of  matrix  systems  (5:33).  The 
structure  of  these  rigid  polymers  is  populated  with  a  high  degree  of  cross-linking  during 
the  curing  process.  As  a  result,  thermosets  cannot  be  reshaped  after  their  original 
formation  and  decompose  thermally  at  high  temperatures  (5:33,  14: 1 1).  The  most 
commonly  used  types  of  thermosets  include  unsaturated  polyesters,  epoxies,  vinylesters, 
and  polyimides  (5:33).  While  epoxies  generally  have  better  mechanical  properties  than 
the  other  thermoset  polymers,  polyimides  are  recognized  for  their  abilities  to  survive 
higher  temperatures. 

Aerospace  Applications 

Applications  of  PMCs  are  often  categorized  based  upon  the  operating  temperature  or 
environment  to  which  they  are  exposed.  PMCs  used  in  aerospace  vehicles  are  likely  to 
be  exposed  to  the  most  demanding  conditions  of  all  applications.  Extreme  temperatures 
encountered  in  these  applications  may  range  from  -150°C  to  550°C  (3:6).  PMCs  may  be 
subjected  to  rapid  heating  conditions,  exposure  to  pollutants  in  the  environment,  solar 
radiation,  and  numerous  other  environmental  factors  (3:6).  PMCs  are  desirable  because 
they  provide  a  lightweight  and  high-strength  solution  to  material  demands. 
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The  need  for  high  temperature  polymer  matrix  resins  to  be  used  in  PMCs  has  been 
evident  for  many  decades.  A  new  class  of  polyimides,  Polymerization  of  Monomeric 
Reactants  (PMR),  was  bom  out  of  the  need  for  a  high  temperature  resistant  polymer 
matrix  resin  that  could  surpass  the  performance  of  the  widely  used  epoxy  resins  in 
aerospace  applications  (20).  Epoxy  resins,  although  frequently  used  in  missile  and 
aircraft  structures,  could  not  sustain  temperatures  above  200°C.  Furthermore,  these  resin 
systems  were  highly  affected  by  environmental  factors  such  as  moisture  which  further 
degraded  their  capabilities.  These  downfalls  led  to  intensive  research  over  the  years  for 
polymers  with  better  elevated  temperature  properties  and  more  environmental  resistance 
than  existing  resins.  In  the  mid-1970’s,  NASA  Lewis  Research  Center  developed  PMR- 
15,  a  neat  matrix  resin  which  offered  stable  performance  at  high  temperatures  in  addition 
to  low  cost  and  ease  in  processing  (16). 

PMR- 15  is  one  of  the  most  sought  after  thennosetting  polyimide  resins  for  high- 
temperature  polymer-matrix  composite  applications  (17:2979).  PMR- 15  polyimide 
composites  have  a  glass  transition  temperature  of  348°C.  Therefore,  the  long-tenn  use 
temperature  of  PMR- 15  composites  is  288°C  (composite  materials  must  be  used  at 
temperatures  well  below  the  Tg  of  the  matrix)  (6).  In  comparison  to  other  resins  of  its 
type,  “PMR- 15  displays  the  best  overall  balance  of  processing,  behavior,  oxidative 
stability,  and  retention  of  mechanical  properties”  at  high  temperatures  (17:2979).  Since 
PMR-15’s  adoption  into  the  commercial  world,  its  processing  materials,  molding 
powders,  chopped-fiber-molding  compounds,  and  variations  of  PMR  resin  are  readily 
available  from  major  aerospace  composite  material  suppliers  (20).  PMR- 15  has  gained 
wide  acceptance  in  the  engineering  community.  It  is  now  being  used  to  fabricate  various 
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engine  components  ranging  from  small  compression  molded  bearings  to  large  structural 
autoclave-molded  engine  ducts  used  on  the  F404  engine  for  the  U.S.  Navy  F-18A  Hornet 
(20).  Additional  examples  of  engine  applications  include  gear  box  covers,  transmission 
housings,  external  exit  flaps,  and  parts  in  numerous  General  Electric  engines  such  as  inlet 
particle  separators,  vent  tubes,  fan  stators,  and  core  cowls  (20).  Today,  PMR-15  has 
been  recognized  as  the  leading  polymer  matrix  resin  for  carbon-fiber-reinforced 
composites  used  in  aircraft  engine  components. 

The  amount  of  polymer  matrix  composites  used  in  fighter  vehicles  has  continuously 
increased  over  the  last  30  years.  For  example,  in  the  1970s  fighter  aircraft  such  as  the  F- 
14,  F-15,  and  F-l  1 1  had  2  to  4  percent  by  weight  carbon-epoxy  composites  (3:1 1).  By 
the  1990s,  these  numbers  grew  to  15  to  30  percent  in  aircraft  such  as  the  A-6,  AV-8,  F- 
18,  and  F-22  by  weight  (3:1 1).  Typical  uses  of  PMCs  have  been  horizontal  and  vertical 
stabilizers,  fuselage  sections,  and  wing  skins.  Similar  trends  are  noticeable  in  rotor 
aircraft.  Whereas  PMCs  were  commonly  used  for  rotor  blades  and  fuel  tanks  in  the 
1970s,  the  V-22’s  entire  airframe  is  composed  of  PMCs  (3:1 1).  The  use  of  lightweight 
polymer  composites  has  increased  the  structural  capability  of  aircraft  while  reducing 
weight,  resulting  in  better  perfonnance.  As  applications  continue  to  mature,  the  use  of 
PMCs  will  likely  be  extended  to  entire  aircraft  and  their  support  structures. 
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Linear  Viscoelastic  Theory 


Viscoelastic  Behavior 

Polymers  exhibit  unique  behavior  under  applied  loads  in  comparison  to  other  matrix 
materials  such  as  metals  or  ceramics.  The  behavior  of  a  material  is  often  characterized  as 
either  an  elastic  solid  or  a  viscous  fluid.  In  the  first  case,  forces  acting  on  a  material 
create  some  deformation  which  is  immediately  reversible  upon  removal  of  the  force  (22). 
Most  materials  which  exhibit  this  behavior  are  linear  elastic  in  that  an  increase  in  stress 
produces  an  equivalent  increase  in  strain  (22).  For  example,  metals  typically  exhibit 
nearly  linear  elastic  behavior  below  their  yield  point.  On  the  contrary,  a  viscous  liquid 
has  no  definite  shape  and  flows  irreversibly  under  the  action  of  external  forces. 

Most  materials  behave  elastically  under  small  stresses  (8:2).  However,  once  the 
stress  exceeds  the  elastic  limit  the  behavior  is  no  longer  elastic  and  the  strain  does  not 
completely  disappear  after  removal  of  the  load  (8:2).  The  permanent  strain  that  remains 
after  the  removal  of  the  load  is  known  as  plastic  strain  (8:2).  Polymers  are  unique  in  that 
they  fit  into  an  intermediate  regime  displaying  the  qualities  of  both  an  elastic  solid  and  a 
viscous  fluid  at  different  times  and  temperatures.  Upon  loading,  they  exhibit  an  initial 
elastic  behavior  followed  by  a  continuous  increase  in  strain  at  a  decreasing  rate.  After 
removal  of  the  load,  they  undergo  an  initial  elastic  recovery  of  strain  followed  by  a 
continuous  decrease  in  strain.  This  type  of  material  response,  which  exhibits  both  solid 
and  liquid  behavior,  is  termed  viscoelasticity  (8:3).  Another  significant  characteristic  of 
viscoelastic  materials  is  dependency  on  time.  Different  loading  rates  used  to  achieve  a 
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certain  stress  level  affect  the  amount  of  strain  accumulated  by  a  material.  Therefore, 
these  materials  are  also  called  time-dependent  materials  (8:2). 

Phenomenological  Aspects 

The  viscoelastic  behavior  of  materials  is  often  studied  with  three  types  of 
experiments  to  include  creep  and  recovery,  relaxation,  and  constant  stress  or  strain  rate 
tests.  Creep  is  defined  as  the  time-dependent  change  in  strain  under  constant  stress 
(28:55).  The  response  to  creep  is  divided  into  three  different  phases:  1)  primary  creep  2) 
secondary  creep  and  3)  tertiary  creep  as  shown  in  Figure  1 .  In  primary  creep,  the  creep 
strain  increases  at  a  decreasing  rate.  During  secondary  creep,  the  creep  strain  rate  is 
nearly  constant.  In  the  final  phase,  tertiary  creep,  the  creep  strain  rate  is  increasing  and 
continues  until  it  terminates  in  fracture.  Generally,  polymers  do  not  enter  into  the  tertiary 
phase  and  only  enter  secondary  creep  if  the  minimum  creep  rate  is  attained  (8:3). 


- Creep  Rate  £ 

Figure  1:  Primary,  Secondary,  and  Tertiary  Creep 
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Following  removal  of  the  load,  recovery  of  strain  occurs  at  a  continually  decreasing 
rate.  For  most  polymers,  a  large  portion  of  the  time-dependent  creep  strain  will  be 
removed  during  recovery  (8:3).  Some  polymers  may  even  exhibit  full  recovery  if 
sufficient  time  is  allowed  (8:3).  Recovery,  like  creep,  is  a  non-linear  process  in  which  the 
rate  is  dependent  on  time.  An  example  creep-recovery  test  is  shown  in  Figure  2  below. 
For  non-linear  viscoelastic  materials,  such  as  polymers,  a  step  input  in  stress  causes  an 
instantaneous  increase  in  strain  followed  by  a  non-linear  accumulation  of  strain 
throughout  the  creep  period.  Similarly,  removal  of  the  load  causes  an  instantaneous 
decrease  in  strain  with  a  continual  non-linear  decrease  in  strain. 


<Jn 

<t0- 


CREEP  RECOVERY 


Figure  2:  Creep-Recovery  Test 

Relaxation  is  defined  as  the  gradual  decrease  in  stress  under  a  constant  strain  as 
shown  in  Figure  3.  The  initial  stress  is  proportional  to  the  applied  strain  and  decreases 
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with  time.  Amorphous  polymers  at  high  temperatures  may  achieve  complete  relaxation 
where  the  stress  eventually  decays  to  zero. 


Figure  3:  Relaxation  Test 


Linearity 

A  material  may  be  considered  linear  viscoelastic  if  it  complies  with  the 
requirements  of  proportionality  and  the  Boltzmann  linear  superposition  principle 
(8:5,  22).  The  proportionality  requirement  is  stated  in  the  following  equation: 

£[ca(t)]  =  c£[a(t)]  (1) 

where  c  is  a  constant,  a  is  the  strain  output,  and  cr  is  the  stress  input.  This  implies  that  a 
stress  input  of  the  scalar  value  c  <j(t )  will  produce  a  strain  output  equivalent  to  that 
produced  by  the  input  <j(t)  multiplied  by  the  constant  c  (8:6,  22).  The  second 
requirement  is  expressed  as  follows: 
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£  [cTj  (t-ti)  +  <J2(t-t2)\  =  S[(JX  ( t  -  t,  )]  +  e[<J2  ( t  - 12 )] 

This  principle  implies  that  if  a  stress,  ax ,  is  applied  at  a  time,  tx ,  and  an  additional 


(2) 


stress,  cr2 ,  is  applied  at  a  later  time,  t2 ,  than  the  strain  output  at  some  time  t  (t  >t2 )  will 
be  equal  to  the  summation  of  the  strain  outputs  from  each  stress  acting  independently 
(8:6,  22).  Figure  4  illustrates  each  of  the  principles  of  linearity.  In  the  superposition 
portion,  Figure  4  (b),  segment  1  denotes  the  strain  due  to  application  of  stress  <jx  ,  segment 

2  denotes  the  strain  due  to  application  of  stress  <j2  ,  and  segment  3  represents  the  strain 
resulting  from  the  summation  of  stress  ax  and  stress  cj2  . 

i)  Proportionality  ii)  Superposition 


a)  b) 

Figure  4:  Principles  of  Linearity 
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The  determination  of  linearity  from  empirical  data  may  be  accomplished  by  the 
formulation  of  isochronous  stress  strain  curves  from  creep  tests  conducted  over  the  stress 
range  of  interest  (22).  The  strain  values  for  different  stress  levels  may  be  selected  at 
various  times  throughout  each  individual  creep  test.  These  values  are  plotted  for  each 
chosen  time  on  a  stress  strain  graph.  If  the  material’s  behavior  is  within  the  linear  region, 
the  curves  will  remain  parallel;  however,  if  the  material  response  becomes  nonlinear  as 
the  stress  values  increase,  the  curves  become  nonparallel.  An  example  of  this  method  is 
shown  in  the  results  section  for  the  PMR-15  neat  resin. 

Mathematical  Representations 

Two  important  material  properties  can  be  defined  for  a  viscoelastic  material 
subjected  to  creep  and  relaxation  tests.  Creep  compliance,  J(t),  is  the  creep  strain  per  unit 
of  applied  stress.  For  linear  viscoelastic  materials,  the  creep  compliance  is  expressed  as 
follows: 

At)  =  ^  (3) 

^0 

The  logarithm  of  creep  compliance  as  a  function  of  the  logarithm  of  time  is  shown 
graphically  in  Figure  5.  During  the  glassy  phase,  J(t)  has  an  order  of  magnitude  of  ~10'9 
Pa'1  (22).  During  the  rubbery  phase,  J(t)  has  an  order  of  magnitude  of  ~10'5  Pa'1  (22). 

The  retardation  time,  or  characteristic  time,  is  the  value  in  the  viscoelastic  region  labeled 
as  the  logarithm  of  f .  This  time  marks  the  inflection  from  rising  to  falling  slope  of  the 
creep  compliance. 
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Figure  5:  Creep  Compliance 


During  creep  testing,  a  polymer  initially  exhibits  “glassy  compliance”  in  which  the  elastic 
deformation  is  a  result  of  bond  deformation  (22).  Eventually  the  compliance  rises  to 
equilibrium  and  exhibits  “rubbery  compliance”  as  a  result  of  molecular  mobility  (8). 

Similarly,  the  relaxation  modulus,  G(t),  is  defined  as  the  relaxation  stress  per  unit  of 
applied  strain.  For  linear  viscoelastic  materials,  the  relaxation  modulus  is  expressed  as 
follows: 

<3(0  =  —  (4) 

^0 

The  logarithm  of  relaxation  modulus  as  a  function  of  the  logarithm  of  time  is  shown 
graphically  in  Figure  6.  During  the  glassy  phase,  G(t)  has  an  order  of  magnitude  of  ~109 
Pa  which  decreases  to  a  value  of  ~105  Pa  during  the  rubbery  phase  (22).  Polymers 
initially  exhibit  “glassy  compliance”  during  a  relaxation  test;  however,  molecular 
mobility  is  the  mechanism  that  generally  prevails  and  is  reached  very  quickly  (22). 
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The  relationship  between  the  creep  compliance  and  relaxation  modulus  for  linear 
viscoelastic  materials  is  expressed  in  the  following  equations: 

G0=-T  G„=J-  (5) 

J  g  Jr 

It  is  important  to  note  that  the  characteristic  time,  or  relaxation  time,  for  relaxation 
response  is  different  than  that  from  the  creep  response  because  relaxation  achieves 
equilibrium  much  faster  (22). 

Integral  Representation  of  Linear  Viscoelasticity 

An  integral  representation  or  hereditary  integral  exists  that  is  capable  of  describing 
creep  strains  due  to  an  arbitrary  stress  input  for  linear  viscoelastic  behavior.  This  integral 
representation  of  viscoelasticity  is  capable  of  incorporating  the  entire  history  of  stress 
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input  into  the  predictive  equation.  This  integral  representation  involves  the  use  of 
Boltzmann’s  Superposition  principle  described  earlier. 

Using  the  unit  step  function  H(t) ,  a  step  input  of  constant  stress  cr0  applied  at  t  =  0 
and  the  corresponding  strain  response  may  be  expressed  as  the  following: 


o-  =  a0H(t) 

(6) 

s(t)  =  a0J(t)H(t ) 

(7) 

If  a  new  stress  input  al  is  applied  not  at  t  =  0  but  at  some  later  time,  c, ,  Eq.  6  and  Eq.  7 
would  become  the  following: 


(8) 


G(t)  =  <TlJ(t-£i)H(t-Z  i) 


(9) 


According  to  Boltzmann’s  superposition  principle,  if  the  stresses  cr0  and  ax  were  applied 
at  their  respective  times,  t(>  and  c, ,  the  resulting  strain  output  becomes  the  following: 


s(t)  =  (J0J ( t)H(t )  +  <Tj  J(t  -  &  )H(t  -  &  ) 


(10) 


If  an  arbitrary  stress  input  is  represented  as  the  sum  of  a  series  of  constant  stress  inputs, 
as  shown  in  Figure  7,  the  stress  input  and  associated  strain  response  may  be  represented 


as: 


o-(t)  =  Ao-;/-/(7-^) 


i= 1 


(ID 

(12) 
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4l  4  2  4  3  £-1  4m  1 

Figure  7:  Arbitrary  Stress  Input 

As  the  number  of  steps  becomes  infinite,  the  summation  may  be  replaced  by  integration 
and  Eq.  12  becomes: 

(13) 

0 

Further  manipulation  to  include  integration  by  parts,  allows  Eq.  13  to  be  written  in  the 
following  form: 

=  (14) 

0 

This  equation  may  be  further  reduced  to: 

m  =  J0a(t)  +  \w-  (15) 

o 
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if  the  compliance  is  separated  into  the  initial  time  independent  elastic  compliance,  J0 , 
and  the  time-dependent  creep  function,  (p(t )  (23:84).  This  same  process  may  be  repeated 
for  a  relaxation  test  in  which  an  arbitrary  strain  is  the  input  and  stress  is  the  response. 

Time-Temperature  Equivalence 

The  concept  of  time-temperature  equivalence  was  introduced  to  allow  linear 
viscoelastic  models  to  provide  more  accurate  predictions.  Empirical  observations  reveal 
that  the  deformation  of  materials,  polymers  in  particular,  often  depends  on  external 
parameters  such  as  temperature  or  prior  aging.  By  plotting  data  on  a  logarithmic  scale, 
however,  it  is  easily  shown  that  often  the  same  process  is  occurring  at  various 
temperatures.  As  a  result,  a  shift  function  is  introduced  to  allow  for  consideration  of 
these  effects.  This  scaling  of  time  to  compare  similar  behavior  is  termed  time- 
temperature  dependence  (22). 

The  time-temperature  dependence  concept  is  demonstrated  in  Figure  8.  In  this 
example,  the  curves  have  identical  shapes  but  differ  by  a  horizontal  shift  factor,  aT  ,  for 
various  temperatures. 
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log(7(/) 


Figure  8:  Dependence  of  Stress  Relaxation  Modulus  on  Time  and  Temperature 


The  time-temperature  dependence  is  described  as  follows: 

l°g(/o)  —  l°g(0  =  l°g(«r)  (16) 

where  t0  is  the  time  corresponding  to  the  master  temperature  T0  and  t  is  the  time 
corresponding  to  some  other  temperature  T .  The  horizontal  shift  is  then  given: 

aT= ^  (17) 

t 

and  the  behavior  at  temperature  T  is  described  as  the  behavior  occurring  at  T0  except  at 
time  t  =  —  .  Because  the  temperature  is  also  a  function  of  time,  this  expression  must  be 
rewritten  as  the  following: 


dtn 


aT(T(t  0)) 


or 


dtn 


aT(T(t  0)) 


(18) 


0  T 
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A  material  that  exhibits  time -temperature  equivalence  is  known  as  a  thermorheollogically 
simple  material  (22).  Williams,  Landel,  and  Ferry  developed  a  commonly  used  equation 
for  detennination  of  the  shift  factor,  known  as  the  WLF  equation: 


log(ar) 


-c.cr-Ty) 

C2+(T-Tref) 


(19) 


where  C,  =  17.4  and  C2  =  51.6.  This  equation  is  often  used  for  amorphous  polymers  with 
Tref  =  Tg  and  T  in  the  nearby  region  of  7),  (2 1:65,  22,  28:109). 


Linear  Viscoelastic  Models 

Rheological  Models 

The  linear  viscoelastic  behavior  of  materials  may  be  represented  by  models  which 
describe  stress-strain  relations  with  time.  These  models  are  composed  of  two  basic 
elements,  the  linear  spring  and  the  linear  viscous  dashpot. 


a)  b) 

Figure  9:  a)  Linear  Spring  and  b)  Linear  Viscous  Dashpot 
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The  linear  spring,  shown  in  Figure  9  a),  results  in  an  instantaneous  strain,  s0 ,  after 
application  of  an  instantaneous  stress,  <j0  ,  represented  by  the  equation: 

cre=Ese  (20) 

where  E  is  the  linear  spring  constant  or  the  Young’s  modulus.  The  behavior  predicted  by 
the  linear  spring  is  shown  in  Figure  10. 

CTA  e 

- 

U  t  tt  t 

Figure  10:  Linear  Spring  Behavior 

The  linear  viscous  dashpot  is  represented  by  the  equation: 

de 

av=r/—  =  r/£e  (21) 

at 

where  77  is  the  coefficient  of  viscosity.  Figure  1 1  depicts  the  behavior  predicted  by  a 
linear  viscous  dashpot.  As  shown  in  Figure  1 1  a),  a  step  input  of  constant  stress  results  in 
a  time  dependent  increase  in  strain.  However,  a  step  input  of  constant  strain,  as  in  Figure 
1 1  b),  results  in  an  infinite  spike  in  stress  followed  by  an  immediate  decrease  to  zero 
stress  at  time  t  =  (0+).  The  stress  then  remains  zero  for  the  remainder  of  time.  The  spike 
in  stress  to  an  infinite  value,  represented  by  the  Dirac  delta  function,  is  physically 
unrealistic. 
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Various  combinations  of  the  linear  spring  and  linear  viscous  dashpot  elements  are 
used  in  rheological  models.  Three  popular  models  include  the  Maxwell  Model,  the 
Kelvin- Voigt  Model,  and  the  Standard  Linear  Solid. 

Maxwell  Model 

The  Maxwell  model  consists  of  a  spring  and  a  dashpot  arranged  in  series  as  shown  in 
Figure  12. 


E  T| 


Figure  12:  Maxwell  Model 
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Because  the  elements  are  arranged  in  series,  the  total  stress  and  total  strain  of  the  model 

are  represented  by  the  following  equations: 

Total  Stress:  a  =  cre  =  av  (22) 

Total  Strain:  e  =  se  +  sv  (23) 

where  the  subscript  e  denotes  the  elastic  spring  element  and  the  subscript  v  denotes  the 
viscous  dashpot  element.  The  total  strain  rate  of  the  model  is  then  represented  as  the 
following: 


Total  Strain  Rate:  e  =  e+  <j„ 


(24) 


By  introducing  the  spring  strain  rate,  obtained  by  taking  the  time  derivate  of  Eq.  20,  and 
the  viscous  damper  strain  rate,  expressed  in  Eq.  21,  into  the  total  strain  rate,  the  total 
strain  rate  for  the  Maxwell  model  becomes  the  following: 


<j  <J 

£  =  £e+£v=~^  +  ^~ 
E  rj 


(25) 


Further  consideration  of  the  stress  equivalence  for  the  model,  described  in  Eq.  22,  allows 
simplification  to  the  following  relation: 


£  = 


&  <j 
— +  — 


E  rj 


(26) 


During  creep,  the  specimen  is  subject  to  a  constant  stress  input  and  as  a  result  a  =  0  . 
By  integrating  with  respect  to  time  and  using  the  initial  condition  a  =  cr0  at  t  =  t0,  the 
strain  response  due  to  creep  as  predicted  by  the  Maxwell  model  is  expressed  as: 


£  = 


—t  +  — 
rj  E 


(27) 
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where  is  the  constant  of  integration  determined  from  initial  conditions.  It  is  now 

possible  to  examine  the  ability  of  the  Maxwell  model  to  predict  material  behavior  for 
creep,  a  step  increase  in  stress,  and  relaxation,  a  step  increase  in  strain.  The  Maxwell 
model  accurately  predicts  an  instantaneous  increase  in  strain  upon  application  of  a  step 
increase  in  stress.  However,  it  predicts  a  linear  increase  in  strain  over  time  which  is  not 
correct.  In  addition,  the  model  only  allows  for  recovery  of  the  elastic  component  of 
strain  and  predicts  a  pennanent  strain  due  to  the  viscous  component  of  the  model.  Creep 
strain  as  a  function  of  time,  as  predicted  by  the  Maxwell  Model,  is  shown  in  Figure  13  a). 


a) 


b) 


Figure  13:  a)  Creep  and  b)  Relaxation  predicted  by  Maxwell  Model 
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In  the  case  of  relaxation,  the  specimen  is  subject  to  a  constant  strain  input  andz  =  0  .  Eq. 

26  becomes  0  =  —  +  — .  Both  sides  of  the  equation  are  multiplied  by  77  and  an 
E  rj 


E 
- 1 

integration  factor  of  e  v  is  used.  By  integrating  with  respect  to  time,  the  prediction  of 


stress  response  during  relaxation  by  the  Maxwell  model  is  represented  by  the  following 


equation: 


o’(0  = 


°oE 

V 


(28) 


where  o 0  is  the  constant  of  integration  detennined  from  initial  conditions.  Relaxation 

stress  as  a  function  of  time,  as  predicted  by  the  Maxwell  Model,  is  shown  in  Figure  13  b). 
The  Maxwell  Model  accurately  predicts  relaxation.  As  depicted  in  Figure  13,  the 
Maxwell  model  allows  for  definition  of  a  relaxation  time  as  follows: 

'*=1  (29) 

E 

During  relaxation  stress  rapidly  decreases  for  times  t  <  tR  with  the  majority  of  relaxation 
occurring  before  tR. 

Kelvin- Voigt  Model 

Similar  to  the  Maxwell  model,  the  Kelvin- Voigt  model  also  consists  of  a  spring  and 
a  dashpot;  however,  they  are  arranged  in  a  parallel  as  shown  in  Figure  14. 
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a 


a 


Figure  14:  Kelvin-Voigt  Model 

The  parallel  configuration  of  the  spring  and  dashpot  elements  in  the  Kelvin-Voigt  model 

result  in  the  following  total  stress  and  total  strain  equations: 

Total  Strain:  s  =  £e  =  sv  (30) 

Total  Stress:  <j  =  <je+<rv  (31) 

where  the  subscript  e  denotes  the  elastic  spring  element  and  the  subscript  v  denotes  the 
viscous  dashpot  element.  The  total  strain  rate  of  the  model  is  then  represented  by  the 
following: 

s  =  se  =  £v  (32) 

The  stress  strain  relations  for  the  Maxwell  model  are  derived  by  considering  the  total 
stress  (Eq.  31),  the  linear  spring  strain  rate  (Eq.  20),  the  viscous  damper  strain  rate,  (Eq. 
21),  and  the  total  strain  rate  (Eq.  32)  as  follows: 

a  =  ae  +  (Tv=E£e  +  rj£v  =  E£  +  rj£  (33) 
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Both  sides  of  Eq.  33  are  divided  by  77  to  provide  the  following  relation. 


0  Es 

—  = - Y£ 

rj  rj 


(34) 


E 
— t 

Integration  is  completed  using  an  integration  factor  of  e  n  .  The  prediction  of  strain 


response  during  creep  by  the  Kelvin  Voigt  model  is  represented  by  the  following 


equation: 


_E_ 

*(0  =  ^r(  \-e^)  (35) 

E 


In  order  to  develop  the  strain  response  during  recovery,  it  is  necessary  to  adapt  an 
alternate  method  for  depicting  creep  using  the  Boltzmann  superposition  principle.  The 
loading  history  may  be  thought  of  as  two  different  steps  in  stress  followed  by  periods  of 
creep  depicted  by  segment  1)  and  2)  in  Figure  15.  Segment  1  illustrates  loading  as  a  step 
increase  in  stress  to  <j0  at  time  t  =  0  followed  by  creep  at  <jg .  Segment  2  depicts 

unloading  of  the  stress  at  time  ti  which  is  equivalent  to  applying  a  stress  (-  cr0 )  at  time 
t  =  ti.  Finally,  recovery  at  zero  stress  is  equivalent  to  creep  at  (-  <jo  )  at  time  >  ti. 
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Figure  15:  Creep  Input  for  Recovery  Prediction  by  Kelvin-Voigt  Model 


Using  Eq.  35,  the  strain  response  for  both  segment  1)  and  segment  2)  can  be  shown  in  Eq. 
36  and  37,  respectively: 


=  - Ml-eV) 
E 


(36) 


E  . 

(J  (.  - (t~h) 

e(t)  =  --±(l-e  *  ) 


(37) 


The  summation  of  Eq.  36  and  Eq.  37  provides  the  strain  response  during  recovery  (t  >  ti) 
by  the  Maxwell  model: 


_E_t 

e(t)  =  -±e  ’  [e*  - 1] 
E 


(38) 
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Figure  16:  a)  Creep  and  b)  Relaxation  predicted  by  Kelvin-Voigt  Model 


The  Kelvin  Voigt  model  accurately  predicts  the  shape  of  the  strain  response  during  creep. 
The  strain  is  increasing  at  an  initial  rate  of  — .  It  continues  to  increase  at  a  decreasing 

rate.  The  bulk  of  the  creep  strain  is  accumulated  before  the  retardation  time  tc  which  is 
defined  as  follows: 


(39) 


In  the  Kelvin-Voigt  model,  stress  is  initially  carried  by  the  viscous  element  and  gradually 
transmitted  to  the  elastic  element  which  eventually  dominates  at  infinite  time.  This 
behavior  is  described  as  delayed  elasticity. 
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In  the  case  of  relaxation,  the  specimen  is  subject  to  a  constant  strain  input  ands  =  0  . 


Equation  34  becomes  —  =  — .  Both  sides  of  the  equation  are  multiplied  by  r/  and  the 
77/7 

Kelvin-Voigt  model  simply  reduces  to  Hooke’s  Law.  Therefore,  because  a  constant 
strain  input  produces  a  constant  stress  input,  the  Kelvin-Voigt  model  is  incapable  of 
predicting  relaxation  behavior. 

Table  1  illustrates  the  abilities  of  the  Maxwell  and  Kelvin  Voigt  models  to  accurately 
predict  the  behavior  associated  with  most  viscoelastic  materials.  Those  characteristics 
which  are  achieved  with  each  model  are  marked  with  an  X.  The  Maxwell  model 
describes  the  time-independent  strain  on  loading  and  unloading  as  well  as  the  relaxation 
behavior.  The  Kelvin-Voigt  model  describes  the  delayed  elasticity  and  recovery  of  the 
creep.  Neither  model  is  independently  capable  of  predicting  the  desired  behavior. 


Table  1:  Ability  of  Maxwell  and  Kelvin-Voigt  Models  to  Represent  Viscoelastic 

Behavior 


Model 

Instantaneous 
increase  in 
strain  on 
loading 

Delayed 

elasticity 

Instantaneous 
decrease  in 
strain 

on  unloading 

Time 

Dependent 

Recovery 

Relaxation 

Maxwell 

X 

X 

X 

Kelvin 

Voigt 

X 

X 

31 


Standard  Linear  Solid 


Numerous  multi-element  models  exist  that  provide  greater  accuracy  in  predicting  the 
behavior  of  a  viscoelastic  material.  One  such  model  is  the  Standard  Liner  Solid  (SLS). 
The  SLS  consists  of  a  spring  element  acting  in  series  with  the  Kelvin  Voigt  model  as 
depicted  in  Figure  17. 


>1 

Figure  17:  Standard  Linear  Solid  (SLS)  Model 


When  detennining  the  governing  equations  for  the  SLS  model,  it  is  useful  to  view  it  as  a 
combination  of  two  elements.  The  first  element  is  the  spring  with  modulus  of  elasticity 
of  Ei.  Its  governing  equation  is  as  follows: 

c ei  —  Exse  j  (40) 

The  second  element  is  the  parallel  combination  of  spring  and  dashpot  referred  to  as  the 
Kelvin  element  represented  by  the  following  equations: 

Spring  equation  <je2  =  E2se2  (41) 

Dashpot  equation  an  =  (42) 
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Total  stress  equation  ak  =  a E2  +  cr 


(43) 


Total  strain  equation  sk  =  sE2  =  s 

The  total  stress  and  total  strain  equations  for  the  entire  model  are  expressed: 


Total  Stress:  cr  =  aEl  =  crk 


Total  Strain:  e  =  sEl  +  sk 


The  total  strain  rate  of  the  model  is  then  represented  by  the  following: 


Total  Strain  Rate:  e  =  em+  sk 


(44) 


(45) 


(46) 


(47) 


Using  the  equations  for  each  element,  Eq.  40-44,  and  substituting  them  into  the  total 
strain  rate  equation  for  the  two  elements  connected  in  series,  Eq.  47,  the  following 
relationship  is  determined: 

.  E2  &  ..  E2  a 

£  +  ^£  =  —  +  (1  +  ^)  — 

V  Ex  Ex  7 


(48) 


crn 


With  an  integration  factor  e  n  and  initial  conditions  of  a  =  a0  and  £  =  —  at  t  =  t0 ,  the 


strain  response  due  to  creep  as  predicted  by  the  SLS  model  is  expressed  as: 


£{t)  =  cr0 


Ex  +E2 

V  E\E2  j 


(49) 


The  recovery  equation  for  the  SLS  model  is  developed  by  applying  the  Boltzmann 
superposition  principle  as  done  for  the  Kelvin  Voigt  model  and  depicted  in  Figure  15.  As 
a  result,  the  strain  response  during  recovery  (t  >ti)  for  the  SLS  model  becomes: 


ar 


HO  =  r>  [e  "  -1] 

En 


(50) 
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Figure  18:  a)  Creep  and  b)  Relaxation  predicted  by  the  Standard  Linear  Solid 

Model 


For  relaxation,  the  input  is  constant  strain  input  and.?  =  0  .  Eq.  48  reduces  to 


—  e  = - h  (1 H — -)  — .  An  integration  factor  of  e  n 

V  Ex  E]  i) 


is  applied  to  this  equation.  By 


integrating  with  respect  to  time,  the  prediction  of  stress  response  during  relaxation  by  the 
SLS  model  is  represented  by  the  following  equation: 


The  SLS  is  capable  of  providing  an  approximate  representation  of  polymer  behavior  in 
the  viscoelastic  range  as  shown  in  Figure  18.  It  is  capable  of  modeling  time-independent 
strain  on  loading  and  unloading,  delayed  elasticity,  time  dependent  recovery,  and 
relaxation.  While  the  Standard  Linear  Solid  model  does  predict  this  general  behavior;  it 
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is  important  to  note  that  it  is  incapable  of  describing  the  behavior  of  many  viscoelastic 
materials  over  a  large  number  of  variables.  Therefore,  more  complex  models  have  been 
proposed  to  provide  greater  accuracy. 


Nonlinear  Viscoelastic  Theory 

The  aforementioned  theories  provide  an  accurate  representation  of  linear  viscoelastic 
material  behavior  if  the  stresses  and  strains  are  small  (24:295).  However,  this  linear 
range  is  often  only  a  small  portion  of  the  total  material  response  before  failure  (24:295). 
Therefore,  further  considerations  are  necessary  to  incorporate  non-linear  behavior  into 
predictive  models.  For  non-linear  materials,  the  creep  compliance,  J(t),  and  the 
relaxation  modulus,  G(t),  are  no  longer  a  function  of  time  alone.  As  shown  in  the 
relations  below,  they  become  functions  of  stress  and  strain,  respectively. 


(52) 

<J 

GU,s)=a{,) 

(53) 

£ 


As  a  result,  the  Boltzmann  superposition  principle  does  not  directly  apply  for  the 
development  of  constitutive  equations  relating  the  stress  strain  behavior  of  materials. 
However,  the  Boltzmann  superposition  principle  is  often  used  as  a  preliminary  basis  for 
developing  relations  that  are  capable  of  predicting  the  behavior  of  non-linear  materials 
(24:296). 

Ledennan  made  two  important  observations  from  the  plots  of  creep  compliance 
(relaxation  modulus)  as  a  function  of  time  obtained  for  certain  textile  fibers.  First,  the 
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creep  compliance  at  various  stress  levels  overlapped  for  small  values  of  time  near  the 
origin  (22).  This  was  an  indication  of  an  instantaneous  elastic  deformation  in  non-linear 
materials.  Secondly,  delayed  creep  and  recovery  was  a  unique  function  of  stress  for  any 
given  load  (22).  This  led  to  the  development  of  the  Modified  Superposition  Principle 
(MSP),  an  extension  of  the  integral  formulation  for  nonlinear  materials.  Leaderman 
proposed  the  following  equation  to  predict  the  strain  behavior  of  non-linear  materials  as  a 
function  of  time,  t,  and  creep  stress,  a  : 

=  (54) 

E  Jr  og 


Note  that  the  Modified  Superposition  Principle  is  very  similar  to  the  integral  fonnulation 
developed  for  linear  viscoelastic  materials  with  the  exception  of  the  function  of  stress  in 
the  second  tenn. 

Pipkin  and  Rogers  extended  the  integral  formulation  for  non-linear  viscoelastic 
materials  by  introducing  a  non-linear  creep  function  for  the  prediction  of  strain  (19:60, 
22).  This  function  is  given  as  the  following: 


C(t,o)  = 


ds(t,a ) 
da 


(55) 


The  prediction  of  strain,  according  to  Pipkin  and  Rogers,  is  then  expressed  as: 


s(t,a)  =  j 


da(r) 

dr 


C[a(r),t  -  t \dz 


(56) 


Pipkin  and  Rogers  used  the  same  approach  for  the  prediction  of  stress  during  a  stress- 
relaxation  test  by  introducing  a  non-linear  relaxation  modulus. 

Schapery  advanced  the  prediction  schemes  for  non-linear  viscoelastic  materials  even 
further  by  incorporating  a  concept  based  on  time -temperature  equivalence  into  the 
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Modified  Superposition  Principle  (24:296).  Schapery’s  model  consists  of  constitutive 
equations  which  include  material  constants  determined  from  empirical  data.  Schapery 
defined  creep  compliance  as  the  summation  of  initial  elastic  compliance  and  secondary 
transient  compliance,  the  first  and  second  term  respectively  in  the  equation  below: 

J(t)  =  ^-  =  D0+AD(t)  (57) 

a 

Accordingly,  the  creep  strain  for  a  case  of  uniaxial  loading  is  expressed  as  follows: 


£(t)  =  g0D0cr{t )  +  gj  J  A D(y/  -  ^')-^[g2cr(r)]dr 


dz 


(58) 


i 

V(0  =  J 


dt' 


a0[o-(O] 


V 


0  “0 


W=j 


dt' 


floW)] 


(59) 


0  u0 


where  g0,gl,g2,a0  are  material  constants  that  depend  on  stress  level  and  y/  is  the 

reduced  time  (18:246,  24:297).  The  reduced  time  defined  by  Schapery  is  based  upon  the 
concept  of  time -temperature  equivalence.  However,  time-temperature  equivalence 
introduces  a  horizontal  shift,  aT ,  dependent  on  temperature  whereas  the  horizontal  shift 
used  in  Schapery’s  model,  aa ,  is  dependent  on  stress.  It  is  convenient  to  recognize  that  if 
g0  =  gj  =  g2  =aa  =  1,  Eq.  58  becomes  the  strain  response  for  a  linear  material  as 

previously  shown  by  the  Boltzmann  Superposition  principle  (24:297). 

Schapery’s  model  is  capable  of  predicting  the  strain  response  when  stress  is  applied 
through  stepwise  loading.  For  example,  in  a  uniaxial  two-step  loading  condition,  the 
applied  stress  is  given  as  follows: 


<7(0  = 


cra,  0  <t<ta 


\ab,  ta<t<  tb 


(60) 
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Therefore,  the  stress  input  may  be  rewritten  as: 

c(t)  =  [H(t)  -  H{t  -  ta  )]aa  +  [H(t  -  ta )  -  Hit  -  tb  ]crh  (6 1 ) 

where  H(t)  is  defined  as  the  unit  step  function.  By  substituting  the  input  stress  of  Eq.  61 
into  Eq.  58,  the  strain  response  due  to  the  first  loading  step,  cra,  is  expressed  as  follows: 


*c(0  = 


g;B(0)+g,"g2"A d„-L 


CTa,  0  <  t  <  ta 


(62) 


While  the  response  to  the  second  loading  step  becomes  the  following: 


*c(0  =  go-D(0)°'i  +gl 


g\^aAD 


f  t  t~t  ^ 

a  a 


\a° 


+  (g2(Tb-g2(Ta)AD 


rt-0 


la  J 


V  a*  J 


,ta<t<th 


(63) 

where  the  superscripts  in  Eq.  62  and  63  represent  the  constants  corresponding  to  the 
stress  level  (18:247).  For  example,  go  =g'o(crfl)  •  These  formulations  may  be  extended 
for  additional  step  inputs  of  constant  stress. 

If  a  creep-recovery  test  is  perfonned  as  shown  in  Figure  2,  the  stress  input  is  reduced 
to: 


cr{t)  =  [Hit)-Hit-tx  )]<r0 

Substituting  this  expression  into  Eq.  58  results  in  the  following  equations  for  strain 
response  during  creep  and  recovery  respectively: 

^c(0  =  goT>(0)cr0  +g!g,AD0—  <r0,  0  <t<tx 

a0 

^(0  =  g2[A£)(^)-A£)(^-^1)K,  t>tx 

where  i//,  =  —  and  y/  =  —  +  (t  -  tx ) 

«o  ao 


(64) 


(65) 

(66) 
(67) 
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These  equations  may  be  further  broken  down  into  the  instantaneous  increase  in  strain 
upon  addition  of  the  load,  <?0 ;  the  final  strain  value  at  the  end  of  the  creep  period,  ex ;  and 

strain  at  the  beginning  of  recovery,  sx ' ,  resulting  from  the  instantaneous  decrease  in 

strain  upon  removal  of  the  load.  These  relations  are  defined  by  extracting  the  appropriate 
tenns  from  Eqs.  65  and  66  as  follows: 


£0  =  g0D(0)cr0 

(68) 

<?,  =  goD(°)°0  +  gig 2^0  ~a0 

(69) 

aa 

t>tl 

(70) 

Figure  19  is  useful  in  gaining  a  visual  understanding  of  these  variables  and  the  behavior 
associated  with  non-linear  viscoelasticity. 


(7n 

ct0- 


tu  / 1 


CREEP  RECOVERY 


Figure  19:  Creep  and  Recovery  Test  Notation  for  Schapery’s  Model 
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Recognizing  that  for  many  engineering  polymers  compliance  can  be  represented  by  a 
power  law,  Schapery  proposed  that  the  transient  component  of  creep  compliance,  the 
second  term  in  Eq.  57,  be  represented  by  a  power  law  as  well: 

A D{¥)  =  C{¥y  (71) 

where  the  constants  C  and  n  are  independent  of  stress  (18:9,  24:32).  Utilizing  this 
relation  the  expressions  for  creep  and  recovery  strain,  Eq.  65  and  Eq.  66,  become  the 
following: 


u(0  = 


g0£>(0)  +  Cglg2 

— LKl  +  ^o)"- 

Sx 


t 

<?0 

\ao  ) 

(72) 

(AaJ'] 

(73) 

where  Asl  is  the  total  amount  of  strain  accumulated  during  the  creep  process  defined  as 
follows: 


t 

As,  —  s,  Sq  =Cg,g2^7a0 

a„ 


(74) 


and  X  is  the  non-dimensional  reduced  time  as  shown  by  the  following: 


(t-t,) 


(75) 


h 

The  uniaxial  constitutive  equations  capable  of  predicting  the  behavior  of  non-linear 
viscoelastic  materials  as  developed  by  Schapery  are  represented  by  Eq.  72  and  73.  It  is 
possible  to  apply  a  characterization  scheme  to  empirical  creep  and  recovery  data  to 
determine  the  parameters  C,  n,  D(0)  and  g0,gl,g2,aeT  as  functions  of  stress  for  a  given 
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material.  Characterization  of  Schapery’s  model  for  the  PMR-15  neat  resin  is  presented 
in  the  results  portion  of  the  paper. 


Previous  Efforts 

Schapery’s  model  is  a  well-established  method  for  predicting  the  behavior  of  non¬ 
linear  viscoelastic  materials.  At  the  onset  of  studies  utilizing  Schapery’s  model  to  predict 
material  behavior  of  this  nature,  Lou  and  Schapery  were  successful  in  characterizing  the 
constitutive  equations  for  a  non-linear  fiber  reinforced  plastic.  The  results  of  this  effort, 
shown  in  Figure  20,  illustrate  the  greatly  increased  capability  to  predict  material  behavior 
when  using  Schapery’s  model  versus  the  Modified  Superposition  Principle  introduced  by 
Leadennan  (13 :22 1 ).  This  research  also  confirmed  the  validity  of  the  power-law 
representation  for  compliance  (13:221). 


Figure  20:  Lou  and  Schapery:  Viscoelastic  Characterization  (13:221) 
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Peretz  and  Weitsman  (18)  conducted  a  similar  experimentation  on  FM-73  adhesive 
in  an  isothermal  environment.  However,  after  characterization  of  Schapery’s  model 
using  creep  and  recovery  data,  the  theory  was  used  to  predict  the  response  under  two-step 
loading  and  hysteresis  response  under  cyclic  constant  stress  rate  tests.  The  two-step 
loading  tests  consisted  of  the  application  of  stress  to  a  value  of  approximately  50%  of  the 
ultimate  tensile  stress,  followed  by  partial  removal  of  the  load  to  produce  different  levels 
of  reduced  stress.  Results  of  the  model  verification  revealed  that  the  model  was  capable 
of  predicting  creep  and  recovery  strain  (18:258).  Peretz  and  Weitsman  emphasized  that 
the  “discrepancy  between  experimental  results  and  the  prediction  of  non-linear  theory  fall 
within  the  scatter  ranges”  of  creep  data  used  to  develop  the  characterization  scheme 
(18:257). 

In  both  works  previously  described,  “mechanical  conditioning”  of  the  specimen  was 
done  to  reduce  the  effects  of  crack  and  damage  growth  prior  to  model  characterization 
(13:220,  18:246).  This  was  done  by  loading  and  unloading  the  specimen  to  a  percentage 
of  the  ultimate  tensile  strength  until  stress-strain  repeatability  was  observed  (13:246). 

The  need  for  this  pre-conditioning  of  the  sample  materials  was  indicative  of  future  work 
that  would  focus  on  the  inclusion  of  a  damage  function  into  Schapery’s  model. 

Smith  and  Weitsman  modeled  the  mechanical  response  of  a  swirl-mat  composite 
composed  of  randomly  oriented  strands  of  E-glass  fibers  embedded  in  a  urethane  matrix 
(26:301).  The  model  predicted  single-step  creep  and  recovery  using  the  process  set  forth 
in  Schapery’s  model  but  also  incorporated  continuum  damage  and  permanent 
deformation  features  which  were  necessary  to  represent  the  behavior  of  the  composite. 
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The  continuum  damage  theory  accounted  for  severe  microcracking  under  increasing 
loads  in  addition  to  overall  reduction  in  stiffness  (26:301).  Model  predictions  were 
compared  to  experimental  results  of  multi-step  creep  and  recovery  tests.  The  above 
process  incorporated  prediction  methodology  for  material  behavior  within  the 
viscoelastic  regime,  damage  due  to  severe  microcracking,  and  permanent  strains  upon 
load  removal.  The  results  revealed  that  it  is  possible  to  maintain  the  creep  expression 
based  upon  power  law,  developed  by  Schapery,  while  introducing  enhancements  for 
additional  behaviors. 

Numerous  other  investigations  have  successfully  captured  the  non-linear  viscoelastic 
behavior  of  materials  using  variants  of  Schapery’s  model.  These  include  efforts  by  Elahi 
and  Weitsman  on  the  mechanical  response  of  P4  Chopped  Glass/Urethane  Resin 
Composite  (7)  and  Jerina  et  al  on  the  characterization  of  random  fiber  composites  (10). 

In  recent  efforts,  Balaconis  characterized  Schapery’s  model  for  BMI-5250-4  neat  resin 
and  concluded  that  the  model  perfonned  well  when  the  stress  levels  applied  were  in  the 
range  used  for  model  characterization  (2).  In  addition,  when  the  prior  stress  rate  was 
comparable  to  that  used  in  model  characterization  experiments,  the  model  was  capable  of 
predicting  stepwise  creep  behavior  (2). 

Since  its  development,  much  research  has  been  done  on  PMR-15  although  most  of  it 
was  focused  on  the  thermo-oxidative  stability  and  effects  of  aging  (17,  27).  As  a  result, 
most  of  these  efforts  do  not  pertain  directly  to  this  investigation.  However,  Marais  and 
Villotreix  did  complete  analysis  and  modeling  of  the  creep  behavior  of  the  PMR-15  neat 
resin  in  the  linear  region.  The  model  employed  in  this  effort  incorporated  the  retardation 
time  concept  with  a  four  element  rheological  model,  known  as  Burgers  model  (15). 
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Burgers  model  connects  Maxwell  and  Kelvin  models  in  series.  The  model  was  capable 
of  predicting  the  creep  compliance  of  PMR-15  for  short  loading  times  for  the  temperature 
range  considered.  The  WLF  equation  was  applied  and  predictions  were  extended  to 
temperatures  beyond  those  used  in  the  model  formulation;  however,  the  predictions 
deviated  significantly  from  experimental  data  (15).  Westberry  conducted  an 
experimental  investigation  of  the  rate-dependent  behavior  and  short-term  creep  of  PMR- 
15  (29).  Results  revealed  that  rate  dependence  was  not  significant  at  elevated 
temperatures  when  conducting  a  simple  tensile  test.  However,  when  examining  the  effect 
of  prior  stress  rate  on  creep  behavior  proceeded  by  uninterrupted  loading,  significant  rate 
dependence  was  observed. 


Thesis  Objective 

The  objective  of  this  thesis  is  to  characterize  a  nonlinear  viscoelastic  model  based  on 
Schapery’s  fonnulation  using  creep  and  recovery  tests  of  PMR-15  neat  resin  at  288  °C. 
The  model  will  be  used  to  predict  the  experimental  results  obtained  under  single  step  test 
as  well  as  multi-step  creep  test  histories.  The  constitutive  model  predictions  will  be 
compared  to  experimental  results  to  reveal  the  ability  of  the  model  to  represent  the 
material  behavior  while  identifying  any  limitations  of  the  model. 

Further  experimental  investigation  will  be  conducted  to  verify  previous  work 
involving  the  rate  dependent  and  history  dependent  behavior  of  PMR-15.  Effect  of 
loading  rate  on  monotonic  stress-strain  behavior  will  be  explored  in  monotonic  tests  at 
several  constant  stress  rates  at  288  °C.  Effect  of  prior  stress  rate  on  creep  behavior  will 
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be  explored  in  creep  tests  preceded  by  uninterrupted  loading  to  a  target  stress  where 
loading  rate  will  be  changed  from  test  to  test.  The  ability  of  the  model  to  take  into 
account  variation  in  stress  rates  and  prior  loading  history  will  be  explored. 
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III.  Material  and  Specimen 


This  section  will  discuss  in  detail  the  material  under  investigation.  The  actual  test 
specimens  will  then  be  discussed,  including  processing  of  the  material,  specimen  shape 
and  tabbing  of  the  specimens. 

Specimen  Development 

Processing 

The  PMR-15  neat  resin  panels  were  provided  by  the  Air  Force  Research 
Laboratory’s  Materials  and  Manufacturing  Directorate.  A  high  purity  PMR-15  imidized 
foam,  commercially  available  from  Cytec  Industries,  Inc,  was  crushed  into  powder  with  a 
hydraulic  press  (29: 10).  Neat  resin  panels  were  compression  molded  from  the  PMR-15 
powder  with  a  heat  press.  Prior  to  machining  the  specimens,  the  C-scans  of  all  panels 
were  examined  and  any  defects  were  noted.  To  discern  the  impact  of  these  defects,  the 
behavior  of  certain  specimens  with  defects  were  compared  to  those  without.  The 
mechanical  testing  of  the  various  specimens  revealed  that  no  recognizable  pattern  was 
observed  between  the  various  qualities  of  specimens  as  reflected  in  the  C-scans. 

Specimen  Geometry 

A  reduced  gage  section  or  dog  bone  shape  was  chosen  for  test  specimens  to  ensure 
gage  section  failures  during  testing.  Test  specimens  were  machined  from  the  PMR-15 
neat  resin  panels  at  the  Air  Force  Institute  of  Technology  Machine  Shop  by  two  different 
methods.  The  first  method  involved  a  water-jet  machine  in  which  a  highly  pressurized 
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mixture  of  water  and  abrasive  cuts  the  specimens.  In  the  second  method,  specimens  were 
cut  using  a  diamond-saw.  While  the  water-jet  process  required  less  time  in  machining 
the  specimen,  the  diamond-grinding  produced  fewer  defects.  While  defects  in  the  grip 
section  would  not  affect  results,  those  occurring  in  the  gage  section  were  unacceptable. 
Specimens  with  machining  defects  in  the  gage  section  were  not  used  in  this  investigation. 

The  thicknesses  of  the  specimens  were  dictated  by  the  panels  from  which  they  were 
cut.  The  average  thickness  of  panels  used  was  3.81  mm  while  the  range  of  thicknesses 
encountered  in  the  specimens  tested  were  from  3.29  to  4.46  mm.  The  specimen  gage 
length  was  17  mm,  fully  enabling  placement  of  the  extensometer  rods  with  a  12.7  mm 
gage  length.  The  specimens  were  stored  at  room  temperature  in  the  laboratory  air.  The 
specimen  geometry  is  shown  in  Figure  2 1 . 


3.81 

1 
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Figure  21:  Test  Specimen  Geometry 


Specimen  Tabbing 

All  PMR-15  neat  resin  test  specimens  were  tabbed  using  a  glass  fabric/epoxy 
material.  The  purpose  of  the  tabs  is  to  evenly  transfer  the  load  from  the  gripping  wedge 
surface  to  the  test  specimen  and  to  protect  the  gripping  surface  of  the  polymer.  The  glass 
fabric/epoxy  material  is  chosen  because  it  is  a  compliant  material  capable  of  providing 
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the  desired  load  transfer.  In  addition,  this  material  absorbs  the  surface  damage  inflicted 
by  the  hydraulic  wedge  grips.  To  apply  the  tabs  to  the  specimen,  the  grip  surface  area  of 
the  specimen  was  cleaned  free  of  debris.  The  tabs  were  then  coated  with  an  M-Bond 
Catalyst.  After  ample  time  for  the  catalyst  to  set  in,  4  to  6  drops  of  M-Bond  200 
Adhesive  were  applied  to  the  tab.  The  tabs  were  placed  onto  the  specimen  and 
continuous  pressure  was  applied  for  at  least  one  minute  to  allow  a  strong  bond  to  form. 

A  tabbed  specimen  is  shown  in  Figure  22. 


Figure  22:  Tabbed  PMR-15  Specimen 
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IV.  Experimental  Setup  and  Testing  Procedures 


This  section  describes  the  equipment  setup  and  associated  testing  procedures  used 
throughout  experimentation. 

Test  Equipment 

The  mechanical  apparatus  necessary  for  testing  included  the  servo  hydraulic  machine 
equipped  with  a  chilled  water  system,  the  extensometer,  the  computer  software,  and  the 
high  temperature  equipment. 

Servo  Hydraulic  Machine 

A  vertically  configured  810  Material  Test  Systems  (MTS)  servo  hydraulic  machine 
was  utilized  for  all  tests.  The  machine  capacity  was  13.3  kN  (3  kips)  whereas  the 
maximum  load  reached  during  all  testing  of  the  PMR-15  resin  was  0.8  kN  (0.2  kips). 

The  hydraulic  machine  employed  MTS  647. 02B  hydraulic  wedge  grips  capable  of 
applying  a  pressure  of  up  to  2 1  MPa  at  each  end  of  the  specimen.  The  grip  pressure 
throughout  testing  of  the  PMR-15  neat  resin  was  set  to  2.76  MPa.  Each  grip  included  a 
pair  of  flat  wedges  coated  with  surf  alloy  to  provide  additional  friction  at  the  gripping 
surface.  The  purpose  of  the  coating  is  to  prevent  the  specimen  from  slipping  throughout 
mechanical  testing.  Additionally,  each  wedge  was  equipped  with  an  inlet  and  outlet  to 
allow  cooling  water  to  pass  through  the  grips  during  testing  to  prevent  overheating  of  the 
grips  region.  The  cooling  water  was  supplied  by  a  NESLAB  model  HX-75  chilled  water 
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system.  The  chiller  regulated  the  temperature  of  the  cooling  water  and  maintained  it 
between  9  and  24°C.  This  system  was  used  during  elevated  temperature  tests.  Finally,  an 
anti-rotation  device  was  fitted  onto  the  machine  actuator  to  maintain  purely  tensile 
loading.  Figure  23  depicts  the  test  setup  with  the  servo  hydraulic  machine,  the  furnace, 
and  the  extensometer  assembly. 


Figure  23:  3  KIP  MTS  Hydraulic  Machine 


Extensometer 

An  MTS  model  632.53E-14  Axial  Extensometer  for  High-Temperature  Testing  was 
used  for  strain  measurement.  The  extensometer  includes  two  3.5  mm  diameter  alumina 
extension  rods  and  has  a  gauge  length  of  12.7  mm  (0.5  in).  Each  extensometer  rod  end 
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has  a  conical  point  which  allows  it  to  maintain  contact  with  a  dimpled  specimen.  For  this 
reason,  two  indentations  were  placed  on  each  specimen  with  the  standard  dimpling  tool 
provided  by  MTS,  shown  in  Figure  24. 


Figure  24:  MTS  Dimpling  Tool 


The  extensometer  tips  are  placed  into  the  dimples  on  the  specimen  to  ensures  that  no 
slipping  of  the  extensometer  occurs  during  testing.  The  use  temperature  for  the 
extensometer  is  up  to  650  °C  without  air  cooling  or  up  to  1200  °C  with  air  cooling.  It  is 
capable  of  measuring  strains  in  the  range  of  -10  %  to  +20%.  The  extensometer,  with  the 
extension  rod  contacting  the  specimen,  is  held  in  place  by  a  spring-loaded  assembly 
which  places  a  3  N  contact  force  per  rod  onto  the  specimen.  The  extensometer  assembly 
also  includes  a  heat  shield  to  protect  the  extensometer  from  the  test  conditions.  The 
extensometer  assembly  is  shown  in  Figure  25. 
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Figure  25:  MTS  High  Temperature  Axial  Extensometer  Assembly 


Computer  Software 

The  MTS  Teststar  II  digital  controller  was  used  for  input  signal  generation  and  data 
acquisition.  The  MTS  Station  Builder  application  was  used  to  create  a  configuration  file 
specific  to  the  performance  of  the  PMR- 1 5  neat  resin.  This  application  allocates 
controller  resources,  specifically  the  input  of  the  appropriate  tuning  parameters  for  a 
material.  The  MTS  Station  Manager  application  was  used  to  perfonn  routine  test 
operations  such  as  controlling  hydraulic  power,  setting  station  interlock  limits,  manually 
controlling  the  force  or  displacement,  monitoring  the  input  and  output  signals  from  the 
controller,  and  running  tests.  As  a  subprogram  of  the  Station  Manager,  the  Multipurpose 
Testware  (MPT)  application  was  used  to  build  test  procedures  by  defining  activities  and 
assigning  triggering  relationships  between  them.  Such  activities  include  loading  the 
specimen,  heating  the  test  specimen,  maintaining  a  particular  force  level,  and  acquiring 
data.  The  construction  of  various  programs  within  the  MPT  software  allowed  for 
automated  testing  and  data  acquisition.  Figure  26  shows  an  example  of  an  MPT  test 
program  used  during  a  creep  test. 
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Figure  26:  MPT  Creep  Test  Procedure 


The  following  signals  were  recorded  during  each  test:  time  (s),  command 
temperature  (°C  ),  displacement  (mm),  force  command  (N),  force  (N),  and  strain  (m/m). 

It  was  important  to  compare  the  force  command  with  the  actual  force  applied  throughout 
each  phase  of  the  test  to  ensure  the  test  was  accurately  executed.  The  frequency  of  data 
acquisition  was  dependent  on  the  test  history  and  the  resultant  deformation  in  the 
specimen.  During  the  heating  phase,  in  which  the  temperature  of  the  test  specimen  was 
elevated  from  23  to  288  °C,  it  was  unnecessary  to  collect  a  large  number  of  data  points  as 
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the  primary  interest  was  the  overall  thermal  expansion  of  the  material  after  reaching  the 
target  temperature.  On  the  contrary,  during  a  loading  phase  when  the  stress  applied  to  the 
specimen  was  increased  from  zero  to  the  target  creep  stress,  the  material  deformation  was 
of  primary  interest  and  the  frequency  of  data  collection  was  increased.  In  this  case,  the 
data  collection  frequency  was  also  dependent  on  the  loading  rate.  For  example,  a  loading 
rate  of  3  MPa/s  required  data  collection  at  a  much  higher  frequency  than  that  of  0.01 
MPa/s  because  the  changes  in  stress  and  strain  were  much  more  gradual  in  the  latter  case. 
It  was  also  important  to  avoid  collecting  data  at  an  excessive  rate  as  this  caused  great 
difficulty  in  transferring  data  to  a  program  for  interpretation.  The  data  acquisition 
frequency  used  for  various  activities  throughout  testing  is  shown  in  Table  2. 


Table  2:  Data  Collection  Frequencies  used  for  Various  Testing  Activities 


Activity 

Data  Collection 
Frequency 
(time  between 
data  points) 

Heat  up  to  Use 
Temperature  of  288  °C 

1  min 

Load/Unload 

Rate  =  3  MPa/s 

0.01  s 

Load/Unload 

Rate  =  1  MPa/s 

0.01  s 

Load/Unload 

Rate  =  0.01  MPa/s 

1  s 

Creep 

1  s 

Recover  for  <  3  h 

1  s 

Recover  for  >  3  h 

3  s 
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High  Temperature  Equipment 


An  MTS  653  Furnace  was  used  for  elevated  temperature  tests.  The  furnace  has  a 
center  split  design  and  an  operating  range  of  100  to  1400  °C.  The  furnace  could  be 
opened  and  shifted  away  from  the  test  area  during  specimen  and  extensometer  setup  and 
closed  around  the  specimen  for  testing.  Each  furnace  half  contained  a  silicon  carbide 
heating  element  surrounded  by  silica  foam.  The  insulation  material  was  carved  to  allow 
for  a  close  fit  around  the  specimen  and  extensometer  rods.  A  19  mm  portion  of  the  test 
specimen  is  encased  by  the  furnace  and  held  at  the  test  temperature.  The  temperature  is 
controlled  by  an  MTS  409.83  Temperature  Controller  Unit.  This  unit  displays  the 
temperature  of  each  furnace  half  and  the  current  input  temperature.  A  feedback  loop 
between  the  controller  and  a  S-type  thennocouple  located  in  the  right  furnace  half 
maintained  the  furnace  temperature  throughout  testing.  Variation  in  temperature  between 
the  right  and  left  halves  of  the  furnaces  was  limited  to  5  °C.  Figure  27  shows  the 
temperature  controller  and  the  right  portion  of  the  furnace  in  place  next  to  the  specimen 
and  the  extensometer. 
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Experimental  Procedures 


Test  Temperature 

Temperature  calibration  was  perfonned  to  detennine  the  input  temperature  needed  to 
produce  the  desired  specimen  temperature  of  288  °C.  Two  K-type  thermocouples  were 
fixed  onto  a  PMR-15  neat  resin  specimen.  The  thermocouples  were  secured  with  Kapton 
tape  and  wire  to  ensure  contact  between  the  thermocouples  and  the  specimen  throughout 
calibration.  The  specimen  was  carefully  situated  in  the  machine  grips  with  the 
thermocouple  wires  exiting  the  oven  through  the  openings  allowed  for  the  extensometer 
rods.  The  thennocouples  were  then  connected  to  an  Omega  Engineering,  Inc.  OMNI- 
CAL-SA-1 10  portable,  two-channel  temperature  sensor  to  monitor  the  specimen 
temperature  in  both  locations.  The  oven  temperature  was  elevated  to  150  °C  at  a  rate  of 
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1.5  °C/min  and  manually  increased  thereafter.  The  manual  increases  were  accompanied 
by  intermittent  dwelling  periods  of  10  min  to  allow  the  specimen  to  soak  at  the  higher 
temperatures.  This  was  necessary  to  ensure  that  the  target  temperature  was  not 
inadvertently  exceeded.  A  set  point  of  253  °C  was  determined,  which  produced  the 
desired  specimen  temperature  of  288  °C.  After  reaching  the  target  temperature,  the  oven 
was  allowed  to  dwell  for  2.5  h  to  ensure  that  the  specimen  temperature  had  thermally 
equilibrated.  During  this  period,  the  temperature  readings  between  the  top  and  bottom 
thennocouples  were  within  ±  4°C  of  the  desired  temperature. 

Test  Procedures 

Numerous  preparatory  steps  were  completed  prior  to  all  testing.  First,  the  Function 
Generator  application  was  used  to  cycle  the  actuator  in  displacement  control  for  at  least 
30  min.  This  was  necessary  to  allow  the  machine  to  properly  warm  up.  Next,  the 
specimen  gage  section  width  and  thickness  were  measured  and  recorded.  These 
dimensions  were  used  to  calculate  the  load  levels  corresponding  to  the  desired  stress 
levels  used  throughout  testing.  The  MPT  test  procedure  was  then  updated  to  reflect  these 
load  values.  The  dimples  were  then  added  to  the  specimen.  The  technique  for  dimpling 
the  specimen  included  marking  the  extensometer  contact  points  on  the  specimen  and 
gently  tapping  the  indentations  into  the  specimen  with  the  dimpling  tool  and  a  hammer. 
Next,  the  grip  pressure  was  set  at  the  appropriate  value  and  the  temperature  controller 
enabled.  The  specimens  were  aligned  in  the  machine  utilizing  the  guides  on  the  rear  of 
the  grip  area.  The  load  cell  was  then  brought  to  the  zero  displacement  level  and  the  top 
portion  of  the  specimen  was  locked  into  the  grip.  After  zeroing  out  the  load  cell  and 
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switching  to  force  control,  the  bottom  portion  of  the  specimen  was  gripped.  The 
extensometer  was  then  mounted  and  the  oven  was  closed  around  the  specimen.  The 
interlocks  were  then  set  to  disable  the  hydraulic  power  if  displacement  and/or  force 
exceeded  expected  limits.  Finally,  the  displacement  and  strain  channels  were  zeroed  out 
before  beginning  the  test.  The  machine  was  kept  in  force  control  at  zero  force  at  the 
onset  of  testing  to  ensure  that  the  specimen  was  allowed  to  freely  expand  during  the  heat 
up  phase.  Finally,  the  test  description,  specimen  label,  procedure  name,  specimen 
dimensions,  and  the  force  values  used  throughout  the  test  were  recorded  prior  to 
beginning  the  test. 


Test  Descriptions 

All  testing  of  the  PMR-15  specimens  was  conducted  at  the  temperature  of  288°C  in 
laboratory  environments.  Therefore,  the  first  step  in  every  test  procedure  was  to  heat  the 
specimen  up  to  the  elevated  temperature.  This  was  achieved  at  a  rate  of  2°C/s.  In 
addition,  each  specimen  was  allowed  to  soak  at  the  target  temperature  for  30  min  prior  to 
testing.  This  was  necessary  to  ensure  that  the  specimen  was  in  thennal  equilibrium  and 
that  all  thennal  expansion  had  ceased. 

Monotonic  Tension  Tests 

Three  monotonic  tensile  tests  to  failure  were  performed  on  PMR-15  specimens.  The 
first  tensile  test  was  perfonned  in  displacement  control  at  a  rate  of  0.025  mm/s  on  a 
specimen  with  no  dimples.  This  test  was  important  because  further  testing  on  dimpled 
specimens  revealed  that  the  fracture  typically  originated  from  one  of  the  dimple 
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locations.  The  extensometer  was  not  used  and  strain  data  was  not  collected  for  this  case 


as  the  dimples  are  necessary  to  maintain  appropriate  extensometer  contact.  The 
subsequent  tensile  tests  were  perfonned  in  stress  control  at  rates  of  1  MPa/s  and  0.01 
MPa/s.  The  main  focus  of  the  tensile  tests  was  to  detennine  an  ultimate  tensile  strength 
from  which  the  stress  levels  in  creep  tests  could  be  determined. 

Creep  and  Recovery  Tests 

Two  creep  tests  were  performed  in  stress  control  at  load  up  rates  of  1  MPa/s  and 
0.01  MPa/s.  The  specimens  were  loaded  to  creep  stress  level  at  the  constant  stress  rate, 
held  at  the  creep  stress  for  a  duration  of  6  h,  unloaded  to  zero  stress  at  the  same  constant 
stress  rate,  and  allowed  to  recover  for  18  h  after  removal  of  the  load.  The  creep  stress 
level,  20  MPa,  corresponded  to  44%  of  the  maximum  tensile  stress  reached  by  the 
undimpled  specimen  in  the  tensile  test  described  above.  These  tests  allowed  for  the 
comparison  ofsload ,  the  strain  accumulated  during  load  up  to  the  creep  stress  level;  e  , 

the  strain  accumulated  during  the  creep  period;  sunload ,  the  strain  removed  during 
unloading;  srec ,  the  amount  of  strain  recovered  at  zero  stress;  and  s  ,  the  pennanent  strain 

remaining  after  recovery  was  saturated.  In  addition,  the  loading  and  unloading  moduli 
were  compared  to  detennine  if  significant  damage  was  inflicted  to  the  material  during  the 
creep  period. 

Monotonic  Load  and  Unload  followed  by  Recovery 

Two  monotonic  load  and  unload  tests  were  perfonned  in  load  control  at  rates  of  1 
MPa/s  and  0.01  MPa/s  followed  by  recovery  at  zero  stress  for  12  h.  The  maximum  stress 
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level  reached  was  20  MPa.  The  loading  and  unloading  moduli,  the  recovery  strain,  and 
the  permanent  strain  remaining  after  recovery  were  compared  for  each  rate. 

Stepwise  Creep  Test 

The  stepwise  creep  test  is  schematically  depicted  in  Figure  28.  This  test  involved 
monotonic  loading  and  unloading  at  constant  stress  rates  with  intermittent  creep  periods 
of  1  hour  duration.  Two  stepwise  creep  tests  were  performed  in  load  control  at  rates  of  1 
MPa/s  and  0.01  MPa/s.  Creep  periods  were  introduced  at  15  and  20  MPa.  After  the  final 
unloading  to  zero  stress,  the  specimen  was  allowed  to  recover  for  12  h.  The  results  of 
this  test  demonstrate  the  effect  of  prior  stress  history  on  creep  and  recovery  behavior. 
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Figure  28:  Stress  as  a  function  of  time  for  stepwise  creep  test 

Model  Characterization  Tests 

Numerous  creep-recovery  tests  were  employed  to  characterize  the  nonlinear 
viscoelastic  behavior  of  PMR-15  using  Schapery’s  model.  The  goal  of  these  tests  was  to 
subject  the  specimen  to  an  instantaneous  increase  in  stress  followed  by  creep  at  various 
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stress  levels.  Similarly,  the  removal  of  the  load  would  occur  instantaneously.  Following 
unloading,  the  specimen  is  allowed  to  recover  at  zero  stress  for  some  duration  of  time.  A 
schematic  of  the  creep  and  recovery  tests  used  in  model  characterization  is  presented  in 
Figure  29. 
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Figure  29:  Stress  as  a  function  of  time  for  creep  and  recovery  tests  used  in  model 

characterization  tests 


The  loading  and  unloading  of  the  applied  stress  cannot  be  accomplished  in  a  step¬ 
wise  manner.  The  stress  rate  of  3  MPa/s  was  used  to  limit  the  loading/unloading  time  to 
less  than  15  s.  The  creep  stress  level  corresponded  to  percentages  of  the  maximum 
tensile  stress  of  45  MPa  reached  by  the  undimpled  specimen  in  the  tensile  test.  The 
percentages  were  ~  1 1, 22,  33,  44,  and  55  percent.  Each  creep  and  recovery  test  was 
conducted  three  times. 

All  testing  is  summarized  in  Table  3.  Once  again,  each  test  was  conducted  at  a 
temperature  of  288  °C. 


61 


Table  3:  Test  Types  and  Test  Specifications  Summary 


Test  Type 

Load/Unload 

Rate 

(MPa/s) 

Stress 

Level 

(MPa) 

%  of  UTS 

Creep 

Time 

(h) 

Recovery 
Time  (h) 

Mono  tonic  Tension 

0.025  mm/s 

— 

— 

— 

— 

Mono  tonic  Tension 

1 

— 

— 

— 

— 

Mono  tonic  Tension 

0.01 

— 

— 

— 

— 

Creep  -  Recovery 

1 

20 

~  44 

6 

18 

Creep  -  Recovery 

0.01 

20 

~  44 

6 

18 

Monotonic  Load/Unload 

1 

20 

~  44 

— 

12 

Monotonic  Load/Unload 

0.01 

20 

~  44 

— 

12 

Stepwise  Creep 

1 

15,20,15 

-33,44,33 

1 

12 

Stepwise  Creep 

0.01 

15,20,15 

-33,44,33 

1 

12 

Model  Characterization* 

3 

5 

-11 

0.5 

1.5 

Model  Characterization* 

3 

10 

-22 

0.5 

5 

Model  Characterization* 

3 

15 

-33 

0.5 

5 

Model  Characterization* 

3 

20 

-44 

0.5 

5 

Model  Characterization* 

3 

25 

-55 

0.5 

5 

*  Each  model  characterization  test  was  triplicated. 
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V.  Results  and  Discussion 


Thermal  Strain 

All  tests  were  conducted  at  a  temperature  of  288  °C  and  the  thermal  stain  value,  sth 

was  recorded  for  each  test  specimen.  The  coefficient  of  thennal  expansion,  a  ,  could 
then  be  determined  from  the  following  relationship: 

sth  =  a- AT  (76) 

where  AT  is  the  temperature  change  in  °C.  The  temperature  change  was  determined  as 
the  difference  between  the  final  specimen  temperature  of  288  °C  and  the  initial  specimen 
temperature  of  23  °C,  which  was  the  nominal  ambient  room  temperature  in  all  tests. 

These  results  are  presented  in  Table  4.  The  average  value  for  the  coefficient  of  thennal 
expansion  was  51.50  x  10"6/°C  and  the  associated  standard  deviation  was  2.82  x  10"6. 

The  average  values  of  coefficient  of  thennal  expansion  for  each  panel  used  throughout 
testing  were  well  within  a  standard  deviation  from  the  total  average.  Specimen  F-l  was 
omitted  from  these  calculations  due  to  the  fact  that  it  was  not  dimpled  and  strain  data  was 
not  collected  from  this  test. 
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Table  4:  Linear  Coefficient  of  Thermal  Expansion 


Specimen 

Number 

Panel 

Number 

Test  Type 

£th  (m/in) 

a 

(10'6/°C) 

F-l 

1 

Mono  tonic  Tension 

N/A* 

N/A* 

F-2 

1 

Mono  tonic  Tension 

0.01430 

53.97 

F-3 

1 

Mono  tonic  Tension 

0.01394 

52.60 

F-4 

1 

Model  Characterization 

0.01550 

58.48 

F-5 

1 

Model  Characterization 

0.01512 

57.06 

F-6 

1 

Model  Characterization 

0.01449 

54.68 

F-7 

1 

Model  Characterization 

0.01351 

50.98 

F-8 

1 

Model  Characterization 

0.01358 

51.23 

F-9 

1 

Model  Characterization 

0.01412 

53.28 

F-10 

1 

Model  Characterization 

0.01446 

54.57 

F-l  1 

1 

Model  Characterization 

0.01348 

50.87 

F-12 

1 

Model  Characterization 

0.01339 

50.52 

F-13 

1 

Model  Characterization 

0.01325 

49.99 

F-14 

1 

Model  Characterization 

0.01317 

49.71 

F-15 

12 

Model  Characterization 

0.01495 

56.42 

F-16 

12 

Model  Characterization 

0.01340 

50.57 

F-17 

12 

Model  Characterization 

0.01396 

52.67 

F-18 

12 

Model  Characterization 

0.01401 

52.85 

F-19 

13 

Model  Characterization 

0.01313 

49.56 

F-20 

1 

Model  Characterization 

0.01336 

50.42 

F-21 

2 

Model  Characterization 

0.01231 

46.44 

F-22 

2 

Model  Characterization 

0.01367 

51.58 

F-23 

2 

Model  Characterization 

0.01323 

49.91 

F-24 

2 

Model  Characterization 

0.01333 

50.32 

F-25 

2 

Monotonic  Load/Unload 

0.01313 

49.53 

F-26 

2 

Monotonic  Load/Unload 

0.01355 

51.13 

F-27 

2 

Creep  -  Recovery 

0.01292 

48.75 

F-28 

2 

Creep  -  Recovery 

0.01307 

49.30 

F-29 

2 

Stepwise  Creep 

0.01245 

47.00 

F-30 

2 

Stepwise  Creep 

0.01191 

44.95 

*  specimen  was  undimpled 


64 


Monotonic  Behavior 


Three  monotonic  tensile  tests  to  failure  were  performed  on  PMR-15  specimens  to 
determine  preliminary  values  of  the  ultimate  tensile  strength  (UTS)  and  modulus  of 
elasticity  at  the  use  temperature.  In  addition,  the  effect  of  stress  rate  on  each  of  these 
properties  was  evaluated.  The  mono  tonic  tensile  test  results  are  summarized  in  Table  5. 


Table  5:  Monotonic  Tensile  Test  Results 


Specimen 

Number 

Panel 

Number 

Loading 

Rate 

UTS 

(MPa) 

E 

(GPa) 

ef 

(%) 

Time  to 
Failure 
(min) 

F-l 

1 

0.025 

mm/s 

48.48 

N/A* 

N/A* 

16 

F-2 

1 

1  MPa/s 

30.49 

1.42 

3.65 

0.53 

F-3 

1 

0.01  MPa/s 

30.38 

1.42 

6.13 

51 

*  specimen  was  undimpled,  no  strain  data  was  recorded 


The  first  tensile  test,  performed  in  displacement  control  at  a  rate  of  0.025  mm/s  on  a 
specimen  without  dimples,  produced  an  UTS  of  48.48  MPa.  The  subsequent  tensile  tests, 
performed  in  load  control  at  rates  of  1  MPa/s  and  0.01  MPa/s,  produced  UTS  values  of 
30.49  MPa  and  30.38  MPa,  respectively.  The  value  of  the  UTS  at  room  temperature  was 
not  measured  in  this  experiment;  however,  the  room-temperature  UTS  value  provided  by 
Cytec  Industries  Inc.  is  38.6  MPa.  It  is  apparent  that  the  introduction  of  dimples  into  the 
virgin  specimens  resulted  in  a  37%  decrease  in  the  ultimate  tensile  strength.  This  is  due 
to  the  fact  that  the  dimples  act  as  crack  starters.  However,  it  is  important  to  recognize 
that  the  UTS  measurement  is  not  the  primary  focus  of  this  research,  while  accurate  strain 
measurements  were  of  greater  concern  and  more  importance.  Therefore,  dimples  were 
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introduced  to  all  remaining  specimens  to  ensure  that  the  extensometer  remained  in 
contact  with  the  material. 

The  modulus  of  elasticity  was  1 .42  GPa,  for  both  tests  perfonned  in  load  control 
with  loading  rates  of  1  MPa/s  and  0.01  MPa/s.  This  modulus  was  determined  by 
calculating  the  slope  of  the  stress-strain  curve  in  the  linear  region.  In  comparison,  the 
modulus  of  elasticity  at  room  temperature  provided  by  Cytec  Industries  was  3.9  GPa. 
These  observations  reveal  that  the  modulus  of  elasticity  decreases  by  64%  as  the 
temperature  increases  from  23  to  288  °C.  Similarly,  the  average  value  of  the  failure 
strain,  sf  ,  was  4.89%  for  the  tensile  tests  conducted  at  288  °C  whereas  the  failure  strain 

at  room  temperature  provided  by  Cytec  Industries  was  1.5%.  The  dependence  of  the 
modulus,  E,  and  failure  strain,  s , ,  on  temperature  is  anticipated  for  a  polymeric  material. 

As  temperature  rises,  polymers’  behavior  changes  from  glass-like  to  rubber-like.  At 
lower  temperatures,  the  stiffness  is  dependent  on  stored  elastic  energy  and  the 
defonnation  is  associated  with  small  molecule  displacements.  The  molecular  chains  have 
an  increased  flexibility  at  higher  temperatures  due  to  increased  free  volume. 

The  stress-strain  curves  obtained  for  the  PMR-15  neat  resin  at  288  °C  are  shown  in 
Figure  30.  Both  curves  display  nearly  linear  elastic  behavior  near  the  origin  and  are 
practically  indistinguishable  in  that  range.  The  approximate  departure  from  linearity 
occurs  at  a  stress  of  approximately  8  MPa  and  strain  of  0.5%.  After  transition  to  the 
inelastic  region,  the  stress-strain  curves  obtained  at  each  loading  rate  become  noticeably 
different.  The  strain  values  at  corresponding  stress  levels  for  tests  perfonned  at  different 
loading  rates  are  compared  to  assess  the  rate  dependence.  Figure  30  reveals  that  the 
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strain  value  corresponding  to  30  MPa  is  3.65%  for  the  stress  rate  of  1  MPa/s  and  6.13% 
for  the  rate  of  0.01  MPa/s.  The  slower  stress  rate  of  0.01  MPa/s  causes  a  68%  increase 
in  a f  compared  to  that  achieved  at  a  rate  of  1  MPa/s.  As  expected,  the  effect  of  stress 

rate  on  Young’s  modulus  was  negligible  as  shown  by  the  overlapping  curves  in  the  initial 
loading  portion  of  Figure  30. 
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Figure  30:  Tensile  stress-strain  curves  for  PMR-15  neat  resin  at  288  °C 
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Effect  of  Prior  Stress  Rate  on  Creep  and  Recovery  Behavior 


Two  creep  and  recovery  tests  were  performed  in  load  control  with  loading/unloading 
magnitudes  of  1  MPa/s  and  0.01  MPa/s.  In  these  tests,  the  stress  was  increased 
monotonically  to  the  creep  stress  level  at  the  aforementioned  rates  and  then  held  constant 
for  a  period  of  6  h  (21,600  s).  The  specimen  was  then  unloaded  to  zero  load  at  the  same 
stress  rate  and  allowed  to  recover  for  a  period  of  18  h  (64,800  s).  These  tests  reveal  the 
effect  of  loading/unloading  stress  rate  on  creep  and  recovery  behavior.  Results  are 
summarized  in  Table  6. 


Table  6:  Effect  of  Prior  Stress  Rate  on  Creep  and  Recovery  Response 


Loading 

Rate 

(MPa/s) 

^ load 

(%) 

f 

creep 

(%) 

^ unload 

(%) 

Srec 

(%) 

(%) 

1 

1.30 

6.67 

1.55 

2.19 

4.20 

0.01 

1.82 

4.28 

1.17 

1.40 

3.53 

As  seen  in  Table  6,  a  decrease  in  the  loading  rate  by  two  orders  of  magnitude  results 
in  a  40%  increase  in  strain  accumulated  during  the  load  up  period,  sIoad .  These  results  are 

congruent  with  those  obtained  in  the  monotonic  tensile  tests  which  revealed  that  a  slower 
monotonic  loading  rate  to  a  desired  stress  level  results  in  a  larger  accumulation  of  strain. 
Next  it  is  important  to  examine  the  influence  of  prior  stress  rate  on  creep  strain 
accumulation,  £c  .  Creep  strain  accumulations  of  6.67  and  4.28%  were  obtained  in 

tests  conducted  with  prior  stress  rates  of  1  MPa/s  and  0.01  MPa/s,  respectively.  Results 
in  Figure  3 1  demonstrate  that  the  creep  behavior  was  considerably  influenced  by  the  prior 
stress  rate.  Westberry  conducted  a  similar  investigation  to  examine  the  influence  of  prior 
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stress  rate  on  creep  strain  at  25  MPa  (29:27).  Creep  strain  accumulations  of  2.83,  2.06, 
and  1.45%  were  obtained  in  tests  conducted  with  prior  stress  rates  of  0.75,  0.075,  and 
0.0075  MPa/s,  respectively.  These  results  support  Westberry’s  qualitative  findings  on 
PMR- 1 5  which  revealed  that  “creep  strain  accumulation  increases  nonlinearly  with 
increasing  prior  stress  rate”  (29:27).  Discrepancies  in  quantitative  results  may  be 
attributed  to  the  fact  that  the  specimens  used  for  testing  in  each  case  were  cut  from  panels 
manufactured  at  different  times.  These  panels  were  likely  formed  from  different  batches 
of  PMR- 15  resin.  In  addition,  the  panels  were  subjected  to  different  post-cure  heat 
treatments  that  may  have  influenced  the  material  properties. 
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Figure  31:  Effect  of  prior  stress  rates  of  1  MPa/s  and  0.01  MPa/s  on  6  hour  creep 

strain  as  a  function  of  time 
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In  order  to  determine  the  creep  phases  that  PMR-15  undergoes,  the  creep  rate  is 
determined  as  a  function  of  time.  Comparison  of  results  in  Figure  32  with  the  creep 
phases  depicted  in  Figure  1  reveals  that  PMR-15  neat  resin  tested  at  288  °C,  with  a  prior 
loading  rate  of  1  MPa/s,  exhibits  primary  and  secondary  creep.  Primary  creep  occurs 
over  the  first  5  h  of  the  creep  test.  This  is  evident  as  the  creep  strain  rate  is  decreasing 
throughout  this  period.  During  the  secondary  creep  phase  the  creep  strain  rate  is  nearly 
constant.  Similar  analysis  of  the  creep  behavior  associated  with  a  prior  loading  rate  of 
0.01  MPa/s  reveals  that  it  transitions  to  secondary  creep  phase  at  approximately  4  h. 
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Figure  32:  Creep  Strain  and  Creep  Rate  as  a  function  of  time  for  6  hour  creep  test 

with  prior  loading  of  1  MPa/s 
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Creep  Rate  (s'1) 


The  secondary  creep  strain  rates  are  calculated  by  examining  a  portion  of  the  creep 
strain  vs.  time  curve  and  plotting  a  best- fit  linear  trend  line  through  the  data.  As  shown 
in  Figure  33,  the  tests  with  different  prior  stress  rates  achieved  secondary  creep  rates  that 
varied  by  an  order  of  magnitude  with  the  higher  prior  stress  rate  resulting  in  larger  creep 
rate.  Westberry  reported  that  creep  tests  conducted  with  varying  prior  stress  rates 
resulted  in  secondary  creep  rates  of  the  same  order  of  magnitude  (29:28). 


Time  (s) 


Figure  33:  Portions  of  creep  strain  as  a  function  of  time  used  to  determine  the 

secondary  creep  rates. 
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After  completion  of  the  creep  tests,  the  specimens  were  unloaded  to  zero  stress  and 
allowed  to  recover  for  12  h.  The  strain  removed  during  unloading,  smload ,  does  not 

indicate  a  strong  dependency  on  the  unloading  rate.  For  the  1  MPa/s  case,  the  amount  of 
strain  removed  during  unloading  was  more  than  the  amount  of  strain  originally 
accumulated  during  loading.  However,  for  the  0.01  MP/s  case,  the  amount  of  strain 
removed  during  unloading  was  less  than  that  gained  during  loading. 

Finally,  it  is  possible  to  examine  srec ,  the  amount  of  strain  recovered  at  zero  stress; 

2lxi&£ p  ,  the  permanent  strain  remaining  after  recovery  for  each  stress  rate.  The  value  of 

£reC  was  measured  as  2.19%  for  the  faster  prior  stress  rate  and  1.40%  for  the  slower  rate; 

however,  it  is  necessary  to  evaluate  these  values  with  consideration  of  the  strain  values  at 
the  initiation  of  the  recovery  period.  For  that  reason,  the  recovery  strain  is  presented  as  a 
percentage  of  initial  strain  at  the  beginning  of  the  recovery  phase.  The  recovery  phase 
begins  after  complete  removal  of  the  applied  stress  during  unloading. 
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Figure  34:  Recovered  strain  shown  as  percentage  of  the  initial  strain  value  at 
beginning  of  recovery  period  following  creep  at  20  MPa  for  various  prior  stress 

rates. 


Figure  34  shows  that  the  specimen  subjected  to  a  higher  prior  stress  rate  recovers  a 
greater  percentage  of  the  strain  measured  immediately  after  unloading.  Both  specimens 
fail  to  achieve  full  recovery  (100%)  and  a  permanent  strain,  s  ,  remains  after  18  h  of 

recovery.  The  actual  recovery  curves  are  presented  in  Figure  35.  It  is  apparent  that 
recovery  has  nearly  saturated  after  18  h,  as  the  recovery  curves  approach  asymptotic 
values.  As  a  result,  allowing  more  time  would  not  result  in  significant  additional 
recovery. 
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Figure  35:  Recovery  at  zero  stress  following  creep  at  20  MPa  for  various  prior  stress 

rates. 


The  loading  and  unloading  moduli  obtained  at  each  rate  were  compared  as  an 
additional  means  for  accounting  for  the  permanent  strain.  As  shown  in  Table  7,  the 
loading  moduli  are  37%  and  59%  greater  than  the  unloading  moduli  for  the  1  MPa/s  and 
0.01  MPA/s  loading  rates,  respectively.  Therefore,  the  6  h  creep  test  at  20  MPa  causes  a 
considerable  reduction  in  stiffness,  indicating  that  damage  had  developed  in  the  material. 
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Table  7:  Creep  Test  Loading  and  Unloading  Moduli 


Loading  Rate 
(MPa/s) 

Loading 

Moduli 

(GPa) 

Unloading 

Moduli 

(GPa) 

1 

1.59 

1.16 

0.01 

1.10 

0.69 

The  results  of  the  creep  and  recovery  tests  conducted  on  PMR-15  at  288°C  indicate 
that  stress  rate  has  a  varying  effect  on  the  strain  accumulation  throughout  each  segment 
of  the  test.  The  amount  of  strain  accumulated  during  loading,  sload ,  was  influenced  by  the 

load  up  rate,  with  greater  strain  resulting  from  the  slower  rate.  The  creep  behavior  at  20 
MPa  depended  on  the  prior  stress  rate.  The  amount  of  strain  recovered  during  unloading 
was  negligibly  influenced  by  the  stress  rate.  The  magnitude  of  the  percent  change  in 
strain  recovered  at  zero  stress  increased  with  increasing  magnitude  of  prior  stress  rate. 


Monotonic  Loading  and  Unloading  followed  by  Recovery 

The  effect  of  stress  rate  on  monotonic  loading,  unloading,  and  recovery  behavior  was 
studied  by  loading  the  specimen  to  20  MPa,  unloading  the  specimen  to  zero  stress,  and 
holding  the  specimen  at  zero  stress  level  for  12  h.  This  test  was  performed  using  stress 
rates  of  1  MPa/s  and  0.01  MPa/s  at  288  °C. 

The  stress-strain  curves  obtained  on  the  loading  and  unloading  paths  for  rates  of  1 
MPa/s  and  0.01  MPa/s  are  compared  in  Figure  36  and  Figure  37.  It  should  be  noted  that 
the  modulus  values  presented  in  these  figures  were  calculated  as  the  slope  of  the  stress- 
strain  curves  within  the  linear  ranges.  For  the  faster  stress  rate,  shown  in  Figure  36,  the 
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stress-strain  curve  remained  nearly  linear  until  approximately  3  MPa.  Past  this  stress 
level  the  slope  of  the  curve  slowly  decreases  until  the  maximum  stress  level  of  20  MPa  is 
reached.  During  unloading,  the  stress-strain  curve  is  initially  non-linear  but  follows  a 
nearly  linear  path  as  the  stress  level  decreases  from  13  MPa  to  zero  stress. 


Strain  (%) 


Figure  36:  Loading  and  unloading  of  PMR-15  neat  resin  at  288  °C  at  a  rate  of  1 

MPa/s 


For  the  slower  prior  stress  rate,  shown  in  Figure  37,  the  stress-strain  curve  displayed 
an  extremely  small  linear  range  up  to  approximately  2  MPa.  Thereafter,  the  increase  in 
stress  resulted  in  a  non-linear  increase  in  strain.  Upon  unloading,  as  the  stress  was 
decreasing  the  strain  continued  to  accumulate.  An  additional  0.5%  of  strain  was 
accumulated  during  unloading  from  20  MPa  to  16  MPa.  After  the  stress  reached 
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approximately  15  MPa,  the  strain  level  began  to  decrease.  A  linear  unloading  path  is  not 
reached  until  very  small  stress  levels  (less  than  4  MPa).  This  curved  unloading  is  a 
distinguishable  feature  observed  at  the  slower  stress  rate. 


Strain  (%) 


Figure  37:  Loading  and  Unloading  of  PMR-15  neat  resin  at  288  °C  at  a  rate  of  0.01 

MPa/s 

While  the  loading  and  unloading  to  20  MPa  is  completed  in  a  time  span  of  42  s  for 
the  faster  rate,  the  slower  rate  requires  4000  s  ( 1 . 1 1  h)  to  complete  this  process.  The 
loading  and  unloading  moduli  for  each  test  are  displayed  in  Table  8.  For  the  1  MPa/s 
rate,  the  loading  modulus  is  approximately  one  percent  greater  than  unloading  modulus. 
The  0.01  MPa/s  stress  rate  exhibits  a  27%  decrease  in  modulus  from  loading  to 
unloading. 
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Table  8:  Monotonic  Loading  and  Unloading  Moduli 


Loading  Rate 
(MPa/s) 

Loading 

Moduli 

(GPa) 

Unloading 

Moduli 

(GPa) 

1 

1.48 

1.46 

0.01 

0.80 

0.63 

The  strain  measured  immediately  at  the  end  of  the  unloading  was  0.17%  for  the 
faster  stress  rate  and  2. 12%  for  the  slower  stress  rate.  The  results  of  the  recovery  tests  are 
shown  in  Figure  38.  Again,  the  recovery  strain  is  presented  as  a  percent  of  the  initial 
strain  measured  immediately  after  unloading.  Due  to  the  small  value  of  that  initial  strain, 
full  recovery  (100%)  is  achieved  within  1  hour  of  the  recovery  period  in  the  test  with 
prior  stress  rate  magnitude  of  1  MPa/s.  The  test  conducted  with  a  slower  prior  stress  rate 
achieves  35%  recovery  after  12  h.  A  permanent  strain  of  1.37%  exists  after  the  12  hour 
recovery.  As  shown  in  Figure  38,  after  the  12  hour  recovery  period,  the  recovery  is 
nearly  saturated  as  the  curve  approaches  a  horizontal  asymptote. 
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Figure  38:  Recovered  strain  shown  as  percentage  of  the  initial  strain  at  beginning  of 
recovery  period  following  monotonic  loading  and  unloading  to  20  MPa  at  various 

stress  rates. 


Stepwise  Creep  Test  and  Recovery 

The  objective  of  the  stepwise  creep  test  was  to  demonstrate  the  effect  of  prior  stress 
history  on  the  mechanical  behavior  of  PMR-15.  The  stepwise  creep  test  subjected  the 
specimen  to  creep  at  two  different  stress  levels  for  periods  of  one  hour.  The  stress  levels 
were  15  and  20  MPa;  creep  periods  were  introduced  during  both  loading  and  unloading  at 
the  15  MPa  level.  Following  the  stepwise  creep  test,  the  specimen  was  allowed  to 
recover  at  zero  stress  for  a  period  of  12  h.  This  test  was  conducted  using  stress  rates  of  1 
MPa/s  and  0.01  MPa/s. 
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Strain  (%) 

Figure  39:  Stress-strain  response  in  a  stepwise  creep  test  with  loading  and  unloading 
stress  rate  magnitude  of  1  MPa/s  and  1  hr  creep  periods  compared  to  uninterrupted 

monotonic  tensile  stress-strain  curve 


The  stress-strain  curve  resulting  from  the  stepwise  creep  test  with  a  loading  rate 
magnitude  of  1  MPa/s  is  presented  in  Figure  39.  The  grey  line  in  this  figure  depicts  the 
stress-strain  curve  obtained  in  the  monotonic  loading  test  conducted  with  the  same  stress 
rate.  The  overlapping  of  the  grey  and  black  curves  demonstrates  that  the  specimen 
behaves  as  expected  during  the  first  loading  portion  of  the  stepwise  test.  The  modulus  of 
elasticity  and  stress-strain  curve  are  consistent  with  the  previously  observed  behavior  of 
PMR-15.  Examination  of  the  stepwise  creep  test  reveals  that  positive  strain 
accumulation,  or  positive  creep,  takes  places  for  the  entirety  of  the  creep  periods  on  the 
loading  path.  However,  during  the  third  creep  period  (i.e  15  MPa  creep  test  on  the 
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unloading  path),  the  amount  of  strain  initially  decreases,  then  increases,  resulting  in 
reversal  of  the  creep  rate.  The  initial  behavior  exhibited  by  PMR-15  is  known  as 
negative  creep.  After  the  negative  creep  was  observed,  creep  rate  reversal  took  place  and 
positive  creep  was  observed. 


Figure  40:  Creep  strain  as  a  function  of  time  for  the  stepwise  creep  test  conducted 
with  a  loading  rate  of  1.0  MPa/s.  L=Loading.  U=Unloading. 


The  creep  strain  produced  during  each  creep  period,  including  negative  creep,  is  also 
presented  as  a  function  of  time  in  Figure  40.  It  is  important  to  note  that  accumulated 
creep  strain  for  each  creep  period  actually  decreases  with  each  additional  step.  While  the 
creep  strain  of  1.73%  was  produced  during  the  first  creep  period  at  15  MPa,  a  smaller 
strain  of  1 .42%  was  accumulated  during  the  next  creep  period  conducted  at  a  higher 
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stress  of  20  MPa.  The  third  creep  period  at  15  MPa  on  the  unloading  path  resulted  in 
only  0.12%  additional  strain.  In  addition,  the  negative  creep  is  apparent  during  the  third 
creep  period.  While  the  amount  of  strain  decrease  during  negative  creep  is  relatively 
small  in  comparison  to  the  other  creep  strains,  the  creep  rate  reversal  behavior  is 
extremely  significant.  The  negative  creep  behavior  of  PMR-15  was  previously  reported 
by  Westberry  (29:3 1).  Westberry  observed  negative  creep  behavior  on  multiple 
unloading  steps  in  which  the  amount  of  negative  strain  increased  as  the  creep  stress 
decreased  (29:3 1).  The  phenomena  of  negative  creep  has  been  observed  in  other 
materials  including  BMI  5250-4  neat  resin  (2:57),  Nylon  66  (12:313),  and  Poly  (ether 
imide)  or  PEI  (12:313).  The  creep  rate  reversal  phenomenon  observed  in  PMR-15  was 
previously  reported  in  some  instances  of  unloading  during  a  stepwise  creep  test  by 
Westberry  (29:41).  In  addition,  Khan  observed  similar  rate  reversal  behavior  at  low 
stress  levels  for  an  amorphous  polymer  known  as  polycarbonate  (11 :81-82). 

The  creep  strain  as  a  function  of  time  for  the  stepwise  creep  test  conducted  with  a 
loading  rate  of  0.01  MPa/s  is  presented  in  Figure  41.  Positive  creep  is  observed  during 
each  creep  period  following  the  slower  loading  and  unloading  rate.  In  this  case,  creep 
strain  accumulation  increases  as  creep  strain  increases  from  15  to  20  MPa.  The  third 
creep  at  15  MPa  during  unloading  produces  the  smallest  creep  strain  accumulation.  The 
first  creep  period  at  15  MPa  produced  creep  strain  of  0.84%,  the  second  creep  period  at 
20  MPa  produced  creep  strain  of  1.36%,  and  the  15  MPa  creep  period  during  unloading 
generated  creep  strain  of  0.33%. 
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Figure  41:  Creep  strain  as  a  function  of  time  for  the  stepwise  creep  test  conducted 
with  a  loading  rate  of  0.01  MPa/s.  L=Loading.  U=Unloading. 


The  creep  strain  accumulated  at  20  MPa  during  the  stepwise  creep  test  was  compared 
to  that  accumulated  during  creep  at  the  same  stress  level  preceded  by  uninterrupted 
loading  at  the  same  loading  rate.  Comparisons  of  the  creep  response  were  made  for  the 
prior  loading  rate  of  1  MPa/s  as  shown  in  Figure  42.  The  creep  behavior  of  the  PMR-15 
neat  resin  is  dependent  on  stress  history.  Strain  accumulated  in  the  creep  test  preceded  by 
uninterrupted  loading  was  approximately  2.5  times  that  accumulated  at  20  MPa  in  the 
stepwise  creep  test.  In  addition,  the  shapes  of  the  creep  curves  are  notably  different  for 
each  loading  case.  The  curve  resulting  from  the  stepwise  creep  test  approaches 
secondary  creep  whereas  the  creep  preceded  by  uninterrupted  loading  clearly  remains  in 
the  primary  creep  regime. 
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Time  (s) 


Figure  42:  Creep  strain  as  a  function  of  time  obtained  at  20  MPa  in  the  stepwise 
creep  test  and  a  creep  test  preceded  by  uninterrupted  loading  at  a  stress  rate  of  1 
MPa/s.  Mechanical  strains  at  the  beginning  of  the  creep  test  are  3.04%  (stepwise) 

and  1.30%  (uninterrupted  loading). 


Recovery  behavior  obtained  in  a  creep  and  recovery  test  was  compared  to  that 
obtained  following  a  stepwise  creep  and  recovery  test.  While  the  creep  and  recovery  test 
produced  greater  creep  strain,  it  also  allowed  for  greater  recovery  at  zero  stress  as  shown 
in  Figure  43.  Dependency  on  prior  history  is  again  evident  in  the  recovery  behavior  of 
PMR-15. 
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Figure  43:  Recovery  at  zero  stress  following  creep  at  20  MPa  in  a  creep  and 
recovery  test  preceded  by  uninterrupted  loading  and  a  stepwise  creep  test. 
Recovered  strain  is  presented  as  a  percentage  of  the  initial  strain  measured 
immediately  after  reaching  zero  stress. 


Viscoelastic  Model 

Numerous  creep-recovery  tests  were  employed  to  characterize  the  nonlinear 
viscoelastic  model  developed  by  Schapery.  After  characterization,  the  ability  of  the 
model  to  predict  material  behavior  was  examined.  The  creep  and  recovery  durations  for 
each  characterization  test  are  presented  in  Table  9.  The  tests  were  repeated  three  times 
for  each  creep  stress.  The  schematic  depicting  the  creep  and  recovery  test  notation  for 
Schapery’s  Model  is  again  shown  in  Figure  44  (presented  again  for  convenience). 
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Table  9:  Model  Characterization  Test  Specifications 


Load/Unload 

Rate 

(MPa/s) 

Stress 

Level 

(MPa) 

%  of  UTS 

Creep 

Time 

(h) 

Recovery 
Time  (h) 

3 

5 

~11 

0.5 

1.5 

3 

10 

~22 

0.5 

5 

3 

15 

~33 

0.5 

5 

3 

20 

~44 

0.5 

5 

3 

25 

~55 

0.5 

5 

*  Each  model  characterization  test  was  triplicated. 
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Figure  44:  Creep  and  Recovery  Test  Notation  for  Schapery’s  Model 
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Linearity  Region 

The  first  step  in  the  model  characterization  was  to  determine  the  region  in  which  the 
behavior  was  linear  viscoelastic.  This  is  necessary  in  order  to  detennine  the  threshold 
stress  level  above  which  the  constants gQ,  g1,g2,aa  are  no  longer  equal  to  one  but  are 
dependent  on  the  stress  level.  The  threshold  stress  was  determined  by  examining 
isochronous  stress  strain  curves  obtained  from  results  of  creep  tests  conducted  over  the 
stress  range  of  5  to  25  MPa  at  increments  of  5  MPa.  As  shown  in  Figure  45,  the  strain 
values  at  various  times  throughout  each  individual  creep  test  were  plotted  for  each  stress 
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level.  The  times  selected  ranged  from  250  s  to  1700  s.  The  isochronous  stress  strain 
curves  of  PMR-15  neat  resin  reveal  linear  behavior  for  stresses  below  10  MPa. 


Strain  {%) 


Figure  45:  Isochronous  stress-strain  curves  of  PMR-15  neat  resin  indicate  region  of 

linear  viscoelastic  behavior 


Panel-to-Panel  Variability 

As  previously  described,  the  goal  of  the  creep-recovery  characterization  tests  was  to 
subject  the  specimen  to  an  instantaneous  increase  in  stress  followed  by  creep  at  various 
stress  levels.  Ideally,  the  removal  of  the  load  would  also  occur  instantaneously.  Because 
it  is  not  possible  to  conduct  the  loading  and  unloading  in  a  stepwise  manner,  a  loading 
and  unloading  rate  of  3  MPa/s  was  used  to  limit  the  loading/unloading  time  to  less  than 
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15  s.  During  the  ramp-up  and  ramp-down  portions  of  the  characterization  tests,  the 
Young’s  modulus  was  calculated  and  compared  for  various  specimens.  It  was  apparent 
that  the  modulus  values  obtained  for  specimens  cut  from  the  same  panel  were 
considerably  different  from  those  obtained  for  specimens  cut  from  different  panels.  The 
average  value  of  modulus  of  elasticity  during  loading  was  1 .65  GPa  with  a  standard 
deviation  of  0. 16  GPa.  The  modulus  values  during  loading  and  unloading  for  various 
specimens  are  presented  in 
Table  10. 


Table  10:  Modulus  During  Loading  and  Unloading  for  Various  Specimens 


Specimen 

Number 

Panel 

Number 

Eloading 

(GPa) 

E/Eavg 

Eunloading 

(GPa) 

F-4 

1 

1.66 

1.0052 

1.62 

F-5 

1 

1.61 

0.9749 

1.72 

F-6 

1 

1.68 

1.0173 

1.68 

F-7 

1 

1.78 

1.0779 

1.68 

F-8 

1 

1.67 

1.0112 

1.42 

F-9 

1 

1.64 

0.9931 

1.71 

F-10 

1 

1.69 

1.0234 

1.72 

F-1 1 

1 

1.76 

1.0657 

1.78 

F-12 

1 

1.82 

1.1021 

1.89 

F-1 3 

1 

1.53 

0.9265 

1.53 

F-1 4 

1 

1.56 

0.9446 

1.53 

F-1 5 

12 

1.29 

0.7811 

N/A 

F-1 6 

12 

1.63 

0.9870 

1.38 

F-1 7 

12 

1.43 

0.8659 

N/A 

F-1 8 

12 

1.61 

0.9749 

1.58 

F-1 9 

13 

1.45 

0.8780 

1.23 

F-20 

1 

1.88 

1.1384 

1.81 

F-21 

2 

2.06 

1.2474 

1.55 

F-22 

2 

1.7 

1.0294 

1.59 

F-23 

2 

1.65 

0.9991 

1.38 

F-24 

2 

1.58 

0.9567 

1.33 
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The  specimens  cut  from  panels  12  and  13  exhibited  significantly  lower  stiffness  than 
those  cut  from  panels  1  and  2.  Note  that  specimens  F-15  and  F-17  failed  during  the  creep 
period  and  therefore  did  not  produce  the  required  data  for  the  model  characterization.  In 
order  to  account  for  panel-to-panel  variability  and  reduce  the  effects  of  data  scatter  in  the 
remaining  specimens,  the  creep  and  recovery  strains  were  adjusted  using  the  following 
equation: 


_  ^ loading 

£adj  ~  F  £ 

havg 


(80) 


Three  creep  and  recovery  tests  were  conducted  at  each  stress  level  of  5,  10,  15,  20  and  25 
MPa.  The  strain  values  recorded  during  the  creep  period  for  each  stress  level  were 
averaged.  The  creep  strain  as  a  function  of  time  for  various  stress  levels  is  presented  in 
Figure  46.  The  heavy  lines  represent  the  average  values  while  the  shaded  areas  represent 
regions  of  data  scatter.  Even  after  adjusting  for  different  modulus  values,  the  data  scatter 
is  apparent.  The  data  scatter  increases  significantly  with  stress  level  and  is  the  most 
apparent  at  the  highest  stress  level  of  25  MPa.  In  addition,  the  scatter  bands  surrounding 
each  stress  level  continuously  expand  during  the  creep  period. 
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Figure  46:  Creep  strains  vs.  time  curves  and  the  associated  data  scatter 


Determination  of  n,  C,  and  D(0) 

The  first  step  in  the  characterization  scheme  is  to  define  a  reference  recovery  strain 
curve  known  as  the  master  curve.  The  curve  is  found  by  comparing  empirical  recovery 

data  to  the  predicted  curve.  In  the  linear  range,  assuming  a  a=  — -  =  1 ,  Eq.  78  reduces  to 

Si 

the  following: 


sr(0  =  [(i  +  A)"-Gl)"] 


(81) 


The  variable  n  becomes  the  only  unknown  in  this  equation.  It  is  determined  by  plotting 
the  recovery  strain,  sr(t),  as  a  function  of  reduced  time,  A  ,  on  a  log-log  plot  for  several 
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trial  values  of  n  .  The  shape  of  the  recovery  strain  curve  is  then  compared  to  empirical 
recovery  data  in  the  linear  region. 

For  the  PMR-15  neat  resin,  the  stress  levels  of  5  and  10  MPa  were  determined  to  be 
in  the  linear  region.  Recovery  duration  following  creep  at  the  stress  level  of  5  MPa  was 
1 .5  h  and  all  specimens  tested  at  this  level  recovered  less  than  70%  of  strain  present  at  the 
onset  of  recovery.  At  the  10  MPa  stress  level,  recovery  duration  was  5  h  and  the 
specimens  recovered  over  80%  recovery  of  the  initial  strain,  as  shown  in  Figure  47. 
Specimen  F-12  recovered  90%  of  initial  strain  and  had  the  lowest  permanent  strain.  For 
this  reason,  the  value  of  n  was  determined  using  the  recovery  data  for  this  specimen. 
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Figure  47:  Recovery  curves  at  zero  stress  following  creep  at  10  MPa.  Recovered 
strain  is  presented  as  a  percentage  of  the  initial  inelastic  strain  measured 
immediately  after  reaching  zero  stress. 


As  shown  in  Figure  48,  the  recovery  strain  at  zero  stress  following  creep  at  10  MPa 
was  plotted  as  a  function  of  reduced  time,  A,  on  a  log-log  plot.  Predictions  using  Eq.  8 1 
for  several  trial  values  of  n  ranging  from  0  to  1  were  plotted  on  the  same  graph  and  the 
predicted  curves  were  compared  to  the  empirical  data.  As  a  result,  n  =0.46  was 
determined  to  provide  the  best  fit  to  experimental  results. 
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Figure  48:  Recovery  at  zero  stress  following  creep  at  5  MPa  plotted  on  a  log-log 
scale  for  determination  of  constant  n. 


After  the  variable  n  was  determined,  two  more  constants  were  determined  using 
empirical  data  from  the  creep  portion  of  the  test.  Recalling  that  g0  =  gl  =  g2  =aa  =1  in 

the  linear  range,  the  examination  of  Eq.  77  for  a  given  creep  stress  level,  cr0 ,  reveals  that 
D(0)  and  C  are  the  only  unknowns: 

£c(t)  =  (D(0)  +  Ctn)cr0  (82) 
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Consideration  of  the  creep  strain  at  two  different  times  during  the  creep  period,  t  =  ta  and 
t  =  tb,  where  0  <  ta  <th  <  t] ,  results  in  the  following  expressions: 


period  and  was  determined  to  be  619.4. 

After  detennining  the  material  constants  n,  C,  and  D(0),  it  is  possible  to  plot  creep 
and  recovery  strains  predicted  by  Eqs.  77  and  78  and  compare  the  predicted  strains  to  the 
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empirical  data  used  to  determine  these  constants.  Once  again,  g0 , g, , g2 ,aa=l  in  the 


linear  range.  The  plots,  shown  in  Figure  49  and  Figure  50,  provided  verification  that  the 
constants  were  correctly  determined. 


Time  (s) 


Figure  49:  Creep  strain  as  a  function  of  time,  experimental  results  used  to  determine 
material  constants  compared  to  predicted  creep  in  the  linear  stress  range 
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Figure  50:  Recovery  strain  as  a  function  of  time,  experimental  results  used  to 
determine  material  constants  compared  to  predicted  recovery  following  creep  in  the 

linear  stress  range. 


Determination  of  stress  dependent  constants 

Next,  the  remaining  material  constants  must  be  determined  for  stress  levels 
beyond  the  linear  range.  The  values  for  gj  and  aa  are  determined  from  the  horizontal 
and  vertical  shifts,  respectively,  needed  to  match  the  master  curve  with  the  corresponding 

experimental  recovery  data  at  each  stress  level.  Figure  5 1  shows  a  vertical  shift,  — - ,  of 

g  i 

the  master  curve  to  match  the  experimental  data.  The  upward  shift  corresponds  to  a  value 


„  Af, 

ot  — L  greater  than  one. 
g  i 


X 


Figure  51:  Vertical  shift  of  master  curve  up  to  experimental  recovery  strain  data 

following  creep  at  15  MPa. 


The  total  amount  of  strain  accumulated  during  the  creep  period,  A sx ,  is  known  from 
experimental  data  and  is  shown  by  the  following  relationship: 


Afj  =e(tl~)-e(t0+) 


(88) 


Therefore,  the  constant  g,  may  be  defined  for  each  stress  level  as  follows: 


g  i  = 


A^i(o-0) 


vertical  _  shift  (cj  Q ) 


(89) 


The  constant  g,  is  then  plotted  as  a  function  of  stress  as  shown  in  Figure  52.  For  values 


of  stress  within  the  linear  range,  gj  is  assumed  to  be  one.  However,  a  trend  line  is 
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applied  to  the  values  in  the  non-linear  range  in  order  to  establish  as  a  function  of 


stress. 


Stress  (MPa) 

Figure  52:  Material  constant,  gb  as  a  function  of  stress 

Similarly,  a  horizontal  shift  of  the  experimental  data  in  Figure  5 1  to  the  left 
corresponds  to  a  value  of  aa<  1.  The  value  of  aa  was  determined  by  calculating  the 

differences  in  time  between  the  vertically  shifted  master  curve  and  the  experimental 
recovery  data  for  a  given  strain.  The  average  was  then  calculated  to  detennine  the  value 
of  aa  at  each  stress  level.  The  constant  aCTis  plotted  as  a  function  of  stress  in  Figure  53. 
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Figure  53:  Material  constant,  aa,  as  a  function  of  stress 

The  process  used  to  determine  the  value  of  g2  is  similar  to  that  used  to  previously 
determine  the  constant  C.  In  the  non-linear  range,  the  difference  of  the  creep  strains  at  ta 
and  tb  becomes: 

(90) 

a0 

which  may  be  solved  for  the  quantity  Z),  defined  below: 

D  _  £A.h)-£c{tg )  _  g,g2C  (91) 

1  s  +  n  jti\  n  ^  ' 

v0(tb  -ta)  a o 
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D]  was  calculated  by  computing  the  average  of  values  of  Dx  when  ta  took  on  all  times 


throughout  the  creep  period  while  th  was  the  time  just  prior  to  the  end  of  the  creep  period. 
After  detennination  of  Dx ,  it  is  now  possible  to  solve  Eq.  91  for  g2  as  follows: 


S  2  = 


Ag  o 

giC 


(92) 


As  the  instantaneous  increase  in  strain  resulting  from  a  step  increase  in  stress,  s0 ,  is 
known  from  experimental  data  the  remaining  unknown  variable,  g0 ,  may  be  found  from 
Eq.  68: 


a  = _ 20 _ 

D(0)a0 

The  material  constants  gn  and  g0  are  plotted  as  functions  of  stress  in 
Figure  54  and  Figure  55,  respectively. 


(93) 
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Figure  55:  Material  constant,  g0,  as  a  function  of  stress 

In  summary,  characterization  of  Schapery’s  model  involved  determination  of  the 
constants  n,  D(0 '),  and  C  which  are  independent  of  stress  level.  It  also  required 
determination  of  aa ,  go,  gj ,  and  g?  which  depend  on  stress.  While  exact  values  for  each 
constant  may  be  detennined  from  empirical  data,  the  extension  of  this  data  to  trends 
provides  the  ability  to  predict  material  behavior  at  stress  levels  other  than  those  used  for 
characterization  of  the  model.  For  this  characterization,  all  values  of  stress  are  in  MPa, 
and  strain  is  returned  in  units  of  microstrain. 
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Model  Verification 


The  validity  of  the  nonlinear  viscoelastic  model  was  tested  against  an  assortment  of 
loading  histories.  The  specimens  used  for  the  model  verification  tests  were  cut  from 
panel  2.  The  tests  included  monotonic  loading  and  unloading  under  constant  stress  rates, 
a  single  step  creep  test  followed  by  recovery,  and  a  stepwise  creep  test  followed  by 
recovery.  Each  of  these  tests  was  carried  out  at  stress  rates  of  1  MPa/s  and  0.01  MPa/s. 
These  tests  were  not  used  to  determine  the  material  constants  and  functions  for  the 
model. 

The  predictions  were  generated  using  a  numerical  integrator  written  to  integrate 
Schapery’s  nonlinear  viscoelastic  model  and  calculate  strain  using  stress  as  a  function  of 
time  as  input  (25).  It  is  useful  to  review  the  basic  constitutive  equations  representing 
Schapery’s  model  for  a  uniaxial  case,  as  previously  discussed: 


e(i)  =  goAMO  +  gi  J -(//') -^-[g2a(r)^T 


0- 


dr 


(94) 


v(t)  =  J 

o 


dt' 


i 

vX*)  =  \ 


dt' 


a0[o-(O] 


o  “o 


(95) 


As  an  is  a  function  of  stress  only,  as  opposed  to  having  additional  time  dependence,  the 


expression  of  ( y/  -  y/  )  was  reduced  to - .  The  elimination  of  this  integral 

aCTO(0) 

improved  the  computational  time  of  the  program.  For  each  recorded  time  step,  the 
integrator  evaluated  s(t)  by  integrating  the  convolution  integral  using  a  Runge  Kutta 
fourth  order  method  for  numerical  integration.  This  allowed  the  integral  portion  of  Eq. 
94  to  become  the  following: 
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(96) 


2  M  (tKO  “  V(T))  ( g2  (<r(r))<r(r)  -  g2  (°(?  ~  d  T))°(T  ~dT)) 

T= 0 

The  equations  representing  y/(t)  were  integrated  prior  to  integrating  the  convolution.  The 
functions  g0 ,  gt,  g2,  and  aa ,  were  characterized  as  piecewise-continuous  functions 

consisting  of  a  region  of  linear  viscoelastic  behavior  below  10  MPa,  where  the  functions 
were  evaluated  as  1,  and  a  nonlinear  region  where  they  became  functions  of  stress.  For 
the  test  conditions  with  a  1  MPa/s  loading  rate  an  integration  time  step  of  0. 1  s  was  used. 
For  the  0.01  MPa/s  rate,  an  integration  time  step  of  1  s  was  used.  The  stress  profiles  were 
dependent  on  the  loading  condition  and  were  input  for  each  individual  case.  The 
integrator  output  the  time,  strain,  and  stress  to  a  comma  separated  variable  file  for  further 
processing  at  each  recorded  time. 

Monotonic  Loading  and  Unloading 

The  first  test  used  for  model  verification  was  monotonic  loading  and  unloading  at  a 
constant  stress  rate  of  1  MPa/s.  The  comparison  of  experimental  results  to  the 
predictions  is  presented  in  Figure  56.  The  model  was  capable  of  predicting  behavior  in 
the  linear  range,  below  10  MPa,  but  predictions  progressively  drifted  from  the 
experimental  data  in  the  nonlinear  range.  A  small  “elbow”  can  be  seen  in  the  predicted 
response  at  the  stress  level  where  the  material  functions  transition  from  linear  to 
nonlinear  behavior  upon  loading  and  unloading.  The  model  over-predicted  the  total 
strain  accumulated  during  loading  as  well  as  the  strain  removed  during  unloading.  The 
model  prediction  does  not  follow  the  same  unloading  path;  however,  the  prediction  does 
run  parallel  to  the  experimental  data  during  unloading  in  the  linear  range.  Note  that  in 
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this  particular  case,  the  model  was  capable  of  predicting  full  recovery  which  was 
congruent  with  the  actual  material  behavior.  This  was  not  surprising  as  only  a  small 
amount  of  strain  was  present  upon  unloading.  It  is  important  to  note  that  the  stress  rate 
used  in  this  test  was  of  the  same  order  of  magnitude  as  that  used  in  model 
characterization  tests. 


%  Strain 


Figure  56:  Loading  and  unloading  at  a  constant  stress  rate  of  1  MPa/s,  experimental 
results  compared  to  nonlinear  viscoelastic  model  predictions 


Comparison  of  the  loading  and  unloading  at  a  slower  loading  rate  of  0.01  MPa/s  is 
shown  in  Figure  57.  In  this  case,  the  model  did  not  predict  the  stress-strain  behavior  of 
PMR-15  accurately.  The  slower  rate  could  not  be  accounted  for  and  the  strain  during 
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loading  is  severely  under  predicted.  The  model  predicts  the  “curved”  unloading 
qualitatively,  but  quantitative  predictions  are  inaccurate.  The  model  predicted  full 
recovery  behavior,  which  was  not  achieved  experimentally.  The  discrepancies  between 
experimental  and  predictive  data  indicate  the  model’s  inability  to  account  for  a  low  stress 
rate  of  0.01  MPa/s. 


%  Strain 


Figure  57:  Loading  and  unloading  at  a  constant  stress  rate  of  0.01  MPa/s, 
experimental  results  compared  to  nonlinear  viscoelastic  model  predictions 


Creep  Followed  by  Recovery 

The  next  test  used  for  model  verification  was  a  creep  test  followed  by  recovery.  The 
predictions  for  the  creep-recovery  tests  are  shown  in  Figure  58  and  Figure  59.  During 
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creep,  the  model  prediction  follows  the  general  shape  of  the  experimental  data  but  the 
model  does  not  accurately  predict  the  creep  strain.  The  predicted  strain  departs  from  the 
experimental  data  after  less  than  one  hour  of  creep.  While  this  creep  test  continued  for  6 
h,  the  tests  used  to  determine  the  model  parameters  were  limited  to  creep  periods  of  0.5  h. 
Extension  of  the  data  scatter  bands,  presented  in  Figure  46,  to  times  greater  than  1800  s 
(0.5  h)  may  indicate  that  the  predictions  of  the  nonlinear  viscoelastic  theory  fall  within 
the  variance  of  the  creep  behavior  of  the  material.  Nonetheless,  the  model  does  not 
provide  an  accurate  quantitative  prediction  for  the  long  term  creep  behavior  shown  here. 


Time  (hr) 

Figure  58:  Creep  at  20  MPa  with  a  prior  stress  rate  of  1  MPa/s,  experimental  results 
compared  to  viscoelastic  model  predictions. 
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The  model  over  predicts  the  recovery  strain  following  creep  at  20  MPa  (see  Figure 
59).  While  the  experimental  results  reveal  that  recovery  is  nearly  saturated  after  12  h,  the 
model  predictions  indicate  that  recovery  is  unsaturated.  The  model  predicts  a  much 
lower  permanent  strain  than  that  resulting  from  the  experimental  results. 


Time  (hr) 

Figure  59:  Recovery  following  creep  at  20  MPa  with  a  prior  stress  rate  of  1  MPa/s, 
experimental  results  compared  to  viscoelastic  model  predictions 


The  model’s  ability  to  account  for  creep  behavior  following  a  much  slower  prior 
stress  rate  was  also  evaluated.  Figure  60  displays  the  experimental  results  of  the  20  MPa 
creep  test  preceded  by  a  stress  rate  of  0.01  MPa/s  together  with  the  model  predictions. 
The  black  dashed  line  represents  the  model  prediction  using  the  0.01  MPa/s  prior  stress 
rate  while  the  grey  dashed  line  represents  the  prediction  with  the  1  MPa/s  prior  stress 
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rate.  During  the  creep  period  the  1  MPa/s  and  0.01  MPa/s  predictions  follow  the  same 
path.  Again,  the  prior  stress  rate  cannot  be  accounted  for  in  the  model.  Therefore, 
predictions  of  material  behaviors  affected  by  prior  stress  rate  are  inaccurate. 


Time  (hr) 


Figure  60:  Creep  at  20  MPa  with  a  prior  stress  rate  of  0.01  MPa/s,  experimental 
results  compared  to  viscoelastic  model  predictions. 


Stepwise  Creep  Test 

A  stepwise  creep  test  was  also  used  for  model  verification.  In  this  case,  the 
predictions  were  compared  to  experimental  data  obtained  with  the  stress  rate  of 
magnitude  1  MPa/s  during  loading  and  unloading.  The  specimen  was  loaded  to  15  MPa, 
loaded  to  20MPa,  unloaded  to  15  MPa,  and  unloaded  to  zero  stress  before  recovering.  A 
creep  period  of  1  h  took  place  at  each  stress  level.  For  this  case,  the  equations  in 
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Schapery’s  model  were  extended  to  include  multiple  steps  as  shown  earlier  in  Eqs.  72-75. 
The  experimental  and  predicted  stress-strain  curves  are  shown  in  Figure  61.  While  the 
amount  of  strain  predicted  during  each  creep  period  appears  to  be  slightly  different  than 
the  actual,  the  model  generally  performs  fairly  well  in  this  case.  A  closer  examination  of 
the  model’s  capabilities  is  presented  in  the  following  graphs. 


Strain  (%) 

Figure  61:  Stress-strain  behavior  of  stepwise  creep  test,  experimental  results 
compared  to  viscoelastic  model  predictions 


Figure  62  displays  strain  as  a  function  of  time  for  creep  periods  at  15  and  20  MPa. 
During  both  creep  tests,  the  model  predictions  are  similar  to  the  experimental  data.  The 
model  slightly  under-predicts  the  creep  strain  value  at  15  MPa  but  provides  an  accurate 
prediction  of  creep  strain  at  20  MPa.  The  effect  of  the  stress  profile  input  on  the  model’s 
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predictive  capabilities  was  evaluated.  While  the  predictions  shown  in  Figure  61  were 
made  by  inputting  the  stress  profde  actually  implemented  during  testing,  predictions 
using  an  instantaneous  stress  step  input  were  also  assessed.  The  dashed  black  line  in 
Figure  62  represents  the  step  input  while  the  solid  grey  line  represents  the  stress  ramp  up 
input.  It  is  seen  that  the  form  of  stress  input  has  no  effect  on  the  model  predictions. 


Time  (s) 


Figure  62:  Creep  behavior  at  15  and  20  MPa  on  the  loading  path  of  the  stepwise 
creep  test,  experimental  results  compared  to  viscoelastic  model  predictions 


The  model  again  perfonned  well  in  predicting  the  strain  as  a  function  of  time  for  the 
15  MPa  creep  test  during  the  unloading  step  (see  Figure  63).  It  is  extremely  important  to 
notice  that  not  only  does  the  model  predict  the  negative  creep  behavior  at  the  beginning 
of  this  creep  period;  it  also  predicts  the  creep  strain  rate  reversal.  The  creep  strain  value 
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at  the  completion  of  this  period  was  also  accurate.  A  magnified  view  of  the  negative 
creep  and  creep  rate  reversal  is  shown  in  Figure  64.  Although  the  model  was  not 
completely  accurate  quantitatively,  qualitatively  it  was  very  powerful.  The  ability  of  the 
model  to  represent  this  aspect  of  the  polymer  behavior  is  apparent. 
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Figure  63:  Creep  behavior  at  15  MPa  on  the  unloading  path  of  the  stepwise  creep 
test,  experimental  results  compared  to  viscoelastic  model  predictions 


Figure  64:  Magnified  view  of  creep  behavior  at  15  MPa  on  the  unloading  path  of  the 
stepwise  creep  test,  experimental  results  compared  to  viscoelastic  model  predictions. 


Finally,  the  predictions  of  recovery  behavior  (see  Figure  65)  were  again 
unsuccessful.  The  model  over  predicts  the  amount  of  strain  that  was  actually  recovered 
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following  the  stepwise  creep  test.  Similar  model  predictions  were  made  of  a  stepwise 
creep  test  with  loading  rate  magnitude  of  0.01  MPa/s.  These  results  are  not  presented  as 
it  has  already  been  shown  that  the  model  cannot  account  for  the  effects  of  rate.  The 
model  predictions  for  the  stepwise  creep  test  conducted  with  the  loading  rate  magnitude 
of  0.01  MPa/s  were  not  accurate.  Again,  the  model  over-predicted  the  recovery  strain. 


Figure  65:  Recovery  following  stepwise  creep  with  a  prior  stress  rate  magnitude  of  1 
MPa/s,  experimental  results  compared  to  viscoelastic  model  predictions 

Schapery’s  nonlinear  viscoelastic  model  was  capable  of  predicting  creep  behavior  of 
PMR-15  neat  resin  at  288  °C  with  varying  degrees  of  accuracy.  The  model  was  capable 
of  accounting  for  prior  loading  history  as  demonstrated  in  the  stepwise  creep  test. 
However,  the  model  could  not  account  for  the  effects  of  prior  stress  rate.  One  significant 
deficiency  of  the  model  was  its  inability  to  accurately  predict  the  recovery  behavior  of 
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PMR-15  in  all  verification  tests.  This  was  likely  due  to  the  fact  that  the  material  did  not 
exhibit  full  recovery  in  the  linear  range.  This  behavior  introduced  a  permanent  strain  that 
is  not  accounted  for  in  the  model.  Because  this  material  does  not  achieve  full  recovery, 
even  in  the  linear  range,  it  is  necessary  to  introduce  permanent  strain  into  the  model 
formulation  in  order  to  accurately  predict  recovery  behavior. 
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VI.  Concluding  Remarks 


Conclusions 

Results  of  the  monotonic  tensile  tests  revealed  that  the  stress  rate  did  not  have  a 
significant  effect  on  the  Young’s  modulus  of  PMR-15  at  288  °C.  However,  the  failure 
strains  were  influenced  by  rate.  As  the  stress  rate  increased  by  two  orders  of  magnitude, 
the  failure  strain  increased  from  3.65  %  to  6.13%.  The  results  of  the  creep  and  recovery 
tests  conducted  on  PMR-15  at  288°C  indicated  that  the  creep  response  was  dependent  on 
the  prior  stress  rate,  with  creep  strain  accumulation  increasing  with  increasing  prior  stress 
rate.  The  magnitude  of  the  percent  change  in  strain  recovered  at  zero  stress  also 
increased  with  increasing  magnitude  of  prior  loading  rate.  Observations  made 
concerning  the  monotonic  loading  and  unloading  of  the  material  revealed  the  effect  of 
stress  rate.  The  stress  rate  had  a  considerable  effect  on  the  stress-strain  behavior,  with 
curved  unloading  becoming  pronounced  at  the  slower  stress  rate. 

The  effect  of  prior  loading  history  on  the  behavior  of  PMR-15  was  studied  by 
subjecting  the  material  to  a  stepwise  creep  test.  For  the  loading  rate  of  1  MPa/s,  the 
amount  of  creep  strain  decreased  as  creep  stress  increased  from  15  to  20  MPa.  Negative 
creep  followed  by  reversal  to  positive  creep  was  observed  in  the  creep  test  at  1 5  MPa 
conducted  on  the  unloading  path.  The  magnitude  of  the  negative  creep  strain  was  smaller 
than  that  of  the  positive  creep  periods  observed  during  loading.  At  the  slower  stress  rate 
of  0.01  MPa/s,  only  positive  creep  was  observed  during  each  creep  period.  The  third 
creep  period  imposed  at  15  MPa  during  unloading  produced  the  smallest  creep  strain. 
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Comparison  of  creep  strains  accumulated  during  a  stepwise  creep  test  to  those 
accumulated  during  creep  preceded  by  uninterrupted  loading  indicate  that  the  prior  stress 
history  affects  the  creep  behavior.  Dependency  on  prior  stress  history  was  also  evident  in 
the  recovery  behavior  following  each  of  these  test  histories. 

The  material  functions  and  constants  in  Schapery’s  model  were  determined  for 
PMR-15  at  288  °C.  Verification  of  the  model  was  accomplished  by  comparing  model 
predictions  to  results  of  various  deformation  experiments.  Schapery’s  nonlinear 
viscoelastic  model  was  capable  of  predicting  the  behavior  of  PMR-15  with  varying 
degrees  of  accuracy.  Predictions  were  most  accurate  when  the  stress  rate  was  of  the  same 
order  of  magnitude  of  that  used  in  the  characterization  tests.  The  predictions  of 
monotonic  loading  and  unloading  behavior  showed  considerable  hysteresis  at  higher 
stress  levels.  These  results  also  produced  over  prediction  of  the  maximum  strain  obtained 
during  the  loading  period.  Comparisons  of  model  predictions  to  creep  data  revealed  that 
the  model  qualitatively  predicted  the  experimental  data  but  did  not  accurately  predict  the 
creep  strain.  The  model  failed  to  provide  an  accurate  quantitative  prediction  of  the  long 
tenn  creep  response.  The  model  performed  fairly  well  in  predicting  the  behavior  of 
PMR-15  during  a  stepwise  creep  test  with  creep  periods  of  one  hour.  The  model 
predicted  negative  creep  and  creep  rate  reversal  in  the  creep  test  conducted  during 
unloading.  Finally,  the  model  was  unable  to  predict  the  recovery  behavior  of  PMR-15. 
This  was  likely  due  to  the  large  permanent  strain  measured  even  in  the  linear  range  of  the 
material. 
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Recommendations 


Future  modeling  efforts  of  the  behavior  of  PMR-15  neat  resin  may  involve 
enhancements  of  Schapery’s  Model  or  an  alternative  constitutive  model.  The  inclusion 
of  a  pennanent  strain  function  into  Schapery’s  model  framework  may  provide  a  more 
realistic  description  of  the  material’s  recovery  behavior  and  result  in  increased  accuracy 
of  the  model  predictions.  Another  model  which  shows  great  potential  for  describing 
polymeric  behavior  is  the  viscoplasticity  theory  based  on  overstress  (VBO).  This  model 
is  particularly  attractive  for  modeling  the  defonnation  behavior  of  PMR-15  as  it  does  not 
require  the  assignment  of  a  linear  range  or  a  purely  elastic  region. 

The  defonnation  of  a  material  is  often  dependent  on  external  parameters  such  as 
temperature,  prior  aging,  or  exposure  to  various  environments.  The  incorporation  of 
additional  variables  of  this  nature  into  constitutive  model  formulations  would  greatly 
extend  predictive  capabilities.  Furthennore,  upon  development  of  a  predictive  model  that 
is  fully  capable  of  describing  the  behavior  of  the  neat  resin  at  hand,  modeling  efforts  must 
be  extended  to  represent  the  behavior  of  composite  materials  to  be  used  in  aerospace 
applications. 
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